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ABSTRACT 

In this paper we study the stellar populations of 356 bright, M r ^ -19, Coma cluster mem- 
bers located in a 2 degree field centred on the cluster core using SDSS DR7 spectroscopy. 
We find - 31% of the sample have significant emission in H/3, [OIII]5007, Ha or [NII]6585, 
due to star-formation or AGN/LINER activity. The remaining portion of the sample we de- 
scribe as passive or quiescent. Using line-ratio diagnostics, we find the fraction of galaxies 
displaying AGN/LINER type emission increases with increasing galaxy luminosity while the 
star-forming fraction decreases. For the quiescent galaxies we find strong correlations be- 
tween absorption line index strength and velocity dispersion (a) for CN2, C4668, Mgb and 
H(3. Employing a planar analysis technique that factors out index correlations with a, we find 
significant cluster-centric radial gradients in H/3, Mgb an d C46 68 for the passive galaxies. We 
use state-of-the-art stellar population models (Schia vonll2007l) and the measured absorption 
line indices to infer the single-stellar-population-equivalent (SSP-equivalent) age and [Fe/H] 
for each galaxy, as well as their abundance patterns in terms of [Mg/Fe], [C/Fe], [N/Fe] and 
[Ca/Fe]. For the passive galaxy subsample we find strong evidence for "archaeological down- 
sizing", with age oc tj - 90±0 06 . This trend is shown to be robust against variations in sample 
selection criteria (morphologically early-type vs spectroscopically quiescent), emission-line 
detection thresholds, index velocity broadening corrections and the specific SSP model em- 
ployed. Weaker positive correlations are obtained between a and all other measured stellar 
population parameters. We recover significant cluster-centric radial stellar population gradi- 
ents for the passive sample in SSP-equivalent age, [Mg/Fe], [C/Fe] and [N/Fe]. These trends 
are in the sense that, at fixed velocity dispersion, passive galaxies on the outskirts of the clus- 
ter are 24% ± 9% younger with lower [Mg/Fe] and [N/Fe] but higher [C/Fe] than those in the 
cluster core. We find no significant increase in cluster-centric radial stellar population gradi- 
ents when fitting to a passive galaxy subset selected to cover the cluster core and South- West 
region, which contains the NGC 4839 subgroup. Thus we conclude that the NGC 4839 in-fall 
region is not unique, at least in terms of the stellar populations of bright galaxies. We speculate 
that the more pronounced cluster-centric radial gradients seen by other recent studies may be 
attributed to the luminosity range spanned by their samples, rather than to limited azimuthal 
coverage of the cluster. Finally, for our passive sample we have found an age-metallicity anti- 
correlation which cannot be accounted for by correlated errors. 

Key words: surveys — galaxies: clusters: individual: Coma — galaxies: evolution — galax- 
ies: stellar content 



1 INTRODUCTION 

Rich galaxy clusters have long been the target of studies attempt- 
ing to address the fundamental questions of galaxy formation and 
evolution. They are ideal laboratories for this task as they harbour 
large numbers of galaxies distributed across a wide range in local 
density and at a common distance. The most direct route to answer 
these questions is to observe clusters at a range of redshifts and 
trace out the assembly history of their galaxy populations. Unfor- 
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tunately this method is observationally inefficient, requiring signif- 
icant amounts of telescope time, and of course often relies on data 
of lower signal-to-noise (S/N). An alternative approach is to con- 
duct detailed analysis of the stellar populations of large samples of 
galaxies in nearby clusters and then use these characteristics as a 
function of their mass, environment or other key properties, to try 
and reconstruct their formation and evolutionary channels. 

The latter of the methods described above, which may be 
termed a type of galaxy archaeology, was pioneered by works 
which identified the environmental dependence of galaxy mor- 
phology in clusters dPressled [1 980). confirmed that red galax- 
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ies, predominately of early-type, follow a well-defined sequence 
in colour-magnitude space while later types occupy a blue cloud 
dAaronson et alJl98ll : lBower et alll992 ; IStrateva et al.l200lh . and 
have analysed the scaling relations seen for early-type galaxies 
jKormendvll 19771; iDressler et al.lll987l ; iDiorgovski & Davisll 19871 ; 
iBenderet all 1992b . 

In a similar vein, the comparison of integrated photometric 
and spectroscopic observations of unresolved galaxies with stel- 
lar population models has become an increasingly valuable tool 
for galaxy archaeology. While it is possible, at least partially, to 
break the well known age-metallicity degeneracy by combining 
appropriate sets of broadband optic al and near-infrared colours 
(Jam es et alj2006l : ICarter et al]2 009). the standard technique typi- 
cally relies on using spectral line indices whic h target a small num- 
ber of info rmatio n-rich absorption features dBurstein et al.|[l984l ; 
iRosdl 1985b . with lFaber etail dl985b_ introducing th e now de facto 
standard Lick indices. In particular IWorthevI d 1994b demonstrated 
how effective pairs of such indices can be at disentangling the ef- 
fects of age and metallicity in unresolved optical spectra. 

Even prior to this seminal work, evidence had been gather- 
ing that nearby early-type galaxies had non-solar a-eleme nt abun- 
dances dPeletied[l989l ; IWorthev et al.lll992l ; lTrageJ[l997h . an ef- 
fect that manifests itself as different metallic line indices im- 
plying different metallicities and ages for a given galaxy. Typ- 
ically this has been found to be more prevalent in giant early- 
type galaxies and is considered evi dence for shorter sta r-formation 
episodes in more massive galaxies dThomas et alj|2005b . This con- 
clusion is driven by the enrichment history of the interstellar 
medium out of which the current generation of stars formed. On 
timescales < 1 Gyr this is dominated by Type II supernova, yield- 
ing more a-elements (normally traced via Mg) than Fe relative 
to the solar neighbourhood, but on longer timescales by Type 
la supernova which redress the balance. With this in mind sig- 
nificant effort has recently been devoted to including non-solar 
alpha element abundance patterns in stellar population models 
dTrager et al]|2000d; [Proctor & Sansorr]|2002l ; IThomas et ai1l200l 
l2004l ; ISchiavorj|2007l ; iDotter et alj|2007l) . 

At ~ 100 Mpc Coma is the nearest rich and dense galaxy clus- 
ter and so the stellar populations of its galaxies have been the sub- 
ject of extensive study for some time. Below we review the key 
findings of a selection of relevant studies on or involving Coma 
members. 

Using multifibre spectroscopic observations of 125 early-type 
galaxies with -20.5 < Ms < -16 located in two fields, one cen- 
tred on the cluster c ore and another ~ 40' to the southwest, 
ICaldwell et al] dl993b identified a number of galaxies with what 
they termed "abnormal" spectra in the latter field. In this case they 
defined abnormal to refer to spectra that show signs of recent star 
formation and/or nuclear activity. They also note that these galax- 
ies are closely associated with an area of enhanced x-ray emission, 
now generally accepted to be indicative of the NGC 4839 group 
merging with th e main bulk of the cluster teriel etal.ll 1992b . 

IfergenserJ ([l 999) analysed the stellar population parameters 
of 115 early- type cluster members with Gunn r < 15, although 
only 71 galaxies in her sample had velocity dispersions and all 
relevant indice s meas ured. Using the stellar population models of 
Vazde kis et alj d 1996b she derived a relatively low median age of 
~ 5 Gyr with significant intrinsic scatter (± 0.18 dex) and found 
[Fe/H] ~ 0.08, again with sizable scatter (± 0.19). She also con- 
cluded that age and [Fe/H] were not correlated with luminosity 
or velocity dispersion. However, strong correlations with these de- 
scriptors were found for [Mg/Fe]. 



The spectroscopic survey of Poggian ti et al. I feooih observed 
278 Coma members (257 with no detectable emission by their cri- 
teria), across the exceptionally wide magnitude range of -20.5 < 
Ms < -14, located in two 32.5 x 50.8 arcmin 2 fields, one tar- 
geted on the cluster core and one on the NGC 4839 subgroup. They 
found that metallicity positively correlates with luminosity but with 
a substantial scatter and that, interestingly, the scatter increases for 
fainter galaxies. In addition they identify an age-metallicity anti- 
correlation in any given luminosity bin, althou gh as they do no t 
take into account correlated errors (see e.g. lKuntschner et al.l200lb 
the strength of this co nclusion is somewh at limited. Further still, 
using the same dataset, ICarter et"al] d2002l . hereafer C02) find evi- 
dence for a cluster-centric radial gradient in the stellar populations 
of their sample galaxies. They interpret this trend as a variation in 
metallicity in the sense that galaxies in the cluster core are more 

metal rich than those in the outer regions. 

Contemporaneously with the Poggianti et al. study, iMoord 

J200lh obtained spectroscopy for 87 bright early-typ e galax ies (Ms 
< -17) in the core of the cluster. Using the [Worthevl d 1994b models 
and a multiple hypothesis testing technique, he found this sample 
to have a sizable metallicity spread, -0.55 ^ [Fe/H] ^ +0.92, but a 
uniform 8 Gyr age of formation with an intrinsic scatter of ± 0.3 
dex (4 to 16 Gyr). He identified no correlation between age and 
velocity dispersion but did find [Fe/H]-<j and [M g/Fe]-q relations. 

More recently, studies s uch a s those of IThomas et al.l d2005b 
andlSanchez-Blazquez et af] d2006b have included Coma members 
in their high density samples when attempting to constrain the en- 
vironmental dependence of galaxy stellar populations. The former 
report that all of the stellar population parameters they measure, 
age, total metallicity and [a/Fe], are correlated with velocity dis- 
persion regardless of local density and that the ages of galaxies in 
low density environments appear s ystematically lower than those in 
high density regions. By contrast ISanchez -Blazgue z et al. (2006) 
find a less well defined picture. Only their low density sample 
has a significant metallicity-cr correlation when using Mgb to es- 
timate metallicity, as opposed to a positive correlation regardless 
of metallicity indicator for their high density subset. Even more 
interestingly they observe age-cr, or so-called "downsizing", rela- 
tions for low density regions and the Virgo cluster but not in their 
Coma galaxies. Finally they find no age-metallicity anticorrelation 
for their high density sample. 

Trager et al. (2008, hereafter T08) conducted a detailed analy- 
sis of 12 early-type Coma cluster members in the magnitude range 
-21.5 < Ms < -16.5 based on Keck/LRIS spectra. They concluded 
that their sample is consistent with having a remarkably young uni- 
form SSP-equivalent age of 5.2 ± 0.2 Gyr and find no indication 
of an age-cr relation, a result which they suggest is supported by 
the majority of previous work on the stellar populations of Coma 
early-type galaxies. 

In the most recent study, at the time of writing, focusing on the 
stellar populations of bright Coma members, Mat kovic et al.l I2009T) 
observed 74 early-type galaxies with -22 ^ Ms ^ -17.5 (~ -20.5 
^ Ms ^ -16) in the cluster core. Their sample spans the velocity 
dispersion range 30 ^ a ^ 260 km s~ , 32 galaxies having a ^ 
1 00 km s _ 1 . Due to large uncertainties in their derived ages they are 
unable to confirm an age-cr relation but find that on average lower a 
galaxies are indeed younger with lower metallicity and [a/Fe] than 
cluster members with higher velocity dispersions. 

Finally, while somewhat fainter in luminosity coverage than 
this work, the most detailed s tudy to date focus ing primarily on 
Co ma dwarf galaxies is that of ISmifh et ail d2008b (hereafter S08) 
and ISmith et al.l d2009ab . They observed 89 red cluster members 
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with -18.5 < Ms < -15.75 split across two fields, one targeting 
the cluster core and one a degree to the south-west of the clus- 
ter centre. They found evidence for a strong cluster-centric radial 
gradient in galaxy age with galaxies in the core typically having 
older ages than systems in the outer region. A similar trend is also 
reported for [Mg/Fe], with stronger Mg enhancement seen in the 
cluster core than the outer region, while no such correlation is ob- 
served for [Fe/H]. 

As such, it is clear that even in this well studied cluster the stel- 
lar population parameters of its host galaxies and how they scale 
with other galaxy properties are still somewhat in debate. In this 
current work we aim to resolve many of the issues discussed above 
and to clarify the effects of velocity dispersion and environment on 
bright Coma cluster members in as homogeneous a fashion as pos- 
sible. To this end we employ spectroscopy from th e Sloan Digital 
Sky Surve£l Data Release 7 dAbaz aiian et al ]2009l) and the up-to- 
date stellar population models o flSchiavorj 020071 . hereafter referred 
to as the Schiavon models), which allow the determination of ele- 
mental abundance patterns in addition to SSP-equivalent age and 
metallicity. 

The paper is outlined as follows. In section 2 we briefly 
overview the SDSS observations and detail our sample selection. 
In section 3 we perform emission line detection and measure our 
galaxies' velocity dispersions and absorption line indices. In sec- 
tion 4 we derive and analyse their stellar population parameters as a 
function of environment and galaxy velocity dispersion. In section 
5 we discuss our findings and in section 6 we review our conclu- 
sions. 
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Figure 1. Colour-magnitude diagram for the 417 Coma members in our 
sample (see text for selection criteria). Red and blue filled circles are non- 
emission line and emission line galaxies (see Section lxTt respectively with 
S/N > 20 per A. Black points are galaxies which did not meet our S/N 
cut or, in the case of detected emission, were not cleaned sufficiently (see 
Section im . The black diamonds highlight galaxies which are significantly 
bluer than the majority of our non-emission sample, determined here by 
them having a (g-r) colour more than 3 times bluer than the rms scatter 
of an unweighted fit to the full non-emission sample at their luminosity. 
The black ringed galaxy, GMP 2946, has a number of flags that indicate 
problems with its Petrosian magnitude (see text for details). 



2 THE SAMPLE 
2.1 SDSS Data 

The SDSS uses a dedicated 2.5 m telescope at Apache Point Obser- 
vatory with a large format mosaic CCD camera for imaging and a 
pair of fibre-fed double spectrographs each with 320 fibres. Both 
spectrographs have a blue and red channel that when combined 
cover ~ 3 800-9100 A at a re solution of R = 1850-2200 (for further 
details see I York et alj|2000|) . The spectra are flux calibrated by the 
SDSS spectral reduction pipeline using 16 spectroscopic standards 
on each plate, colour selected to be F8 subdwarfs. The SDSS fibre 
diameter of 3" equates to 1.4 kpc at Co ma assuming Hp = 7 1 km 
s" 1 Mpc~\ Q M = 0.27 and fi A = 0.73 fainshaw et ai]|2009l) . We 
adopt a distance modulus m — M — 35.0 for Coma. 

All spectroscopy for this work comes from the SDSS Main 
Galaxy Sample, th e target selection of which is detailed in 
IStrauss et alj fc002h . Briefly, SDSS imaging is used initially for 
star-galaxy separation and then to measure an r-band Petrosian 
magnitude and define a circular aperture containing half a galaxy's 
Petrosian flux. The Main Galaxy Sample then consists of galaxies 
with r pctr o ^ 17.77 and mean surface-brightness < fi r >so ^ 24.5 
mag arcsec~ 2 with rflbrc < 19, after correcting for Galactic ex- 
tinction. Objects with r fibre < 15 are rejected to limit the effects of 
cross talk when extracting the spectra of neighbouring faint objects, 
as are targets that are flagged as saturated, bright or blended and not 
deblended by the photometric pipeline. For reference, applying the 
above constraints to the SDSS photometric database we estimate ~ 
96 % of galaxies in our Coma field (see below) have SDSS spectra. 
It is more than likely the remaining fraction result from having a 



neighbour within 55", the minimum separation permitted between 
fibres. 



2.2 Sample Selection 

The rationale behind our sample selection is simple, we aim to in- 
clude as many Coma memb ers a s possible and exten d the recent 
work of ISmith et al] d2008l) and ISmith et all {2£>09a) to brighter 
magnitudes. To this end we include all galaxies with r pct ro *S 16 
within 2° (~ 3.3 Mpc) of our designated cluster centre at RA = 
12:59:48, Dec = +27:58:50, which is approximately half way be- 
tween NGC 4889 and NGC 4874, that have spectroscopy in SDSS 
DR7. We note that all SDSS photometric values quoted through- 
out the rest of this paper are Petrosian magnitudes unless otherwise 
specified and have been corrected for Galactic extinction and k- 
correction using the provided parameters from the SDSS database. 
In addition we opt for an inclusive redshift cut of 0.01 ^ z ^ 
0.04 based on ~ 4cr li mits of the cluster velocity dispersion from 
IColless & Dunn] dl99rj) and taking Coma to be at z = 0.0230 
These criteria result in our sample containing 417 confirmed Coma 
members. We note that a more conservative 3a cut only removes 
six galaxies and does not significantly affect our results. 

Our magnitude cut is set by a trade off between maximising 
sample size, and completeness, and obtaining galaxy spectra that 
have sufficient S/N; in this case we opt for a median S/N per A 
20 between 4000-6000 A, to suitably constrain stellar population 
parameters (see Section |4~TV Effectively this amounts to a cut of 
fflbrc i$ 18, which results in ~ 87% (363) of our sample having the 
desired S/N. Of course the majority of our incompleteness, due to 



http://www.sdss.org 



http://nedwww.ipac.caltech.edu/ 



4 J. Price et al. 



Table 1. Classification data for the four blue passive galaxies highlighted 
by black diamonds in Fig. [TJ Morphology is taken from NED and the 
quoted H<5 equivalent widths, in units of A, are those measured by the SDSS 
pipeline. 



ID 


(g-r)o 


Morphology 


EW(H<5) sdss 


GMP 2640 


0.33 


SOp 


7.9 


GMP 3892 


0.44 


SBO/a 


5.6 


GMP 3439 


0.50 


SB0 


2.3 


GMP 4974 


0.51 


SAO 


3.75 



our required S/N, occurs in the final 15.75 r 16.0 bin where 
we are ~ 44% complete. 

In Fig. Q] we plot the colour-magnitude diagram for our sam- 
ple. As detailed in Section [37X1 galaxies are, at least initially, only 
subdivided further based on whether they show detectable emission 
and as such it is interesting to see from Fig. Q] that almost all spec- 
troscopically quiescent galaxies (98%) that meet our selection re- 
quirements are compatible with being on the cluster's red sequence. 
The remaining five galaxies are worth further comment. 

The galaxy highlight ed by a black circle in Fig. Q] is GMP 
2946 ( IGodwin et alj|l983l) . This object is projected onto the halo 
of the giant galaxy NGC 4889, and its SDSS photometry is erro- 
neous. More careful treatment shows that it lies on the cluster's red 
sequence. We tabulate the relevant details of the four other blue pas- 
sive galaxies, denoted by black diamonds in Fig.Q] in TableQ]where 
the morphological information is taken from NED and the H<5 ab- 
sorption line equivalent widths are th ose measured by the SDSS 
pipeline SPECTROlD. We note that the Poggianti et al. (2004]) clas- 
sification requirements for post-starburst k+a galaxies are no emis- 
sion and EW(H<5) > 3 A. Thus three out of four of these galaxies 
are compatible with being k+a systems. Furthermore, all three are 
located toward the faint end of our sample's magnitude range, a 
fact which is perhaps not surprising given the lack of brig ht blue 
k+a galaxies in Coma as reported by Poggi anti et alj J2004). Here- 
after we will refer to these four galaxies as the blue subset of the 
passive sample and where comparisons to literature studies involv- 
ing colour selected samples are necessary, we intend to mask this 
subset. 



3 SPECTROSCOPIC MEASUREMENTS 
3.1 Emission line detection 

To conduct our stellar po pulation analysis w e must first identify 
and reject, or correct (e.g. lTrager et al.|[2000al) . galaxies with neb- 
ular emission which acts to infill the stellar H/3 absorption, our 
key B aimer line age indicator, and so results in overestimated stel- 
lar ages. We could, of course, opt to use the higher order B aimer 
lines H7 or H<5 which in turn suffer less from emission infill. How- 
ever, here we prefer H/3 for its relative stability against abundance 
pattern variations and, hence, its essential role in the fitting code 
we shall employ (see Section l4~lT l. We also wish to avoid includ- 
ing galaxies with emission in [NIJ5198 w hich acts to perturb the 
Mgb index ( Goudfrooii & Emsellern| [l996). Fortunately the broad 
wavelength coverage of the SDSS spectra allows access to Ha and 
[NII]6583 which are ideal for flagging galaxies in our sample that 
host active galactic nuclei (A GN or LINER) or a re currently under- 
going star formation (see e.g. lSmith et al.ll2007h . 

As a first pass we make use of the emission line measurements 
taken by SPECTROlD. The code uses a median/mean filter to fit 
each spectrum's continuum with emission and absorption features 



masked. Emission, and pseudo-Lick absorption, lines are then mea- 
sured on the continuum subtracted spectrum as constrained single 
gaussians, with lines in close proximity to one another being fit as a 
blend. Here we deem galaxies with equivalent widths (EW) in any 
of H/3, [OIII]5007, Ha or [NII]6585 > 0, detected at the 3a level, 
to be emission-line candidates. This cut selects 186 galaxies with 
S/N > 20. 

The method of detecting and measuring emission lines has 
been sig nificantly improved recently with the introduction of GAN- 
DALF bvlSarzi et alj d2006t). itself an extension of the pixel fitting 
code of Cappellari & Emsellem|j2004l) which will be detailed later. 
Instead of requiring the observed spectra to be continuum sub- 
tracted prior to analysis, they show that it is possible to simulta- 
neously fit the stellar component and nebular emission lines. This 
approach acts to remove the inherent bias caused by the masking 
of emission line regions, which often also contain important ab- 
sorption features, when fitting the stellar continuum in the former 
method. Indeed we note that as the technique used by the SDSS 
does not make any attempt to disentangle absorption and emission 
components, we seek to check the robustness of the survey's emis- 
sion line detections. 

As input, GANDALF requires a set of template spectra with 
higher spectral resolution than the observations. Here we use 75 
randomly se lected spectra of ty pes F to M from the Indo-U.S. Stel- 
lar Library ( Valdes et alj 2004), which have a wavelength cover- 
age of 3460-9464 A observed at FWHM ~ 1 A. We also include 
another 10 templates from the same library which are specifically 
chosen to have strong Mgb indices and so further help minimise 
template mismatch in systems with strong magnesium absorption. 
The code then finds the best linear combination of the template 
spectra and emission lines, again parameterised as gaussians, in a 
X 2 sense with respect to the observed data. 

We opt to fit our sample in the wavelength range 4000-6800 
A as this includes Ha and [NII]6585, our primary emission indica- 
tors, and, if detected, allows for the emission features which affect 
the Lick indices relevant to this work to be corrected. Following 
Sarzi et al. I d2006l) we require that bona-fide emission lines be de- 
tected at an amplitude-to-noise ratio (A/N) 4, where noise is de- 
fined as the residual-noise level between the best fit and observed 
spectrum and therefore this cut actively screens against template 
mismatch generating fake emission lines. As before, a galaxy with 
this level of emission in any of H/3, [OIII]5007, Ha or [NII]6585 
is flagged. By comparing EW as a function of A/N for each line, 
we empirically estimate that A/N = 4 corresponds to EW(Ha) ~ 
0.5 A, EW([NII]6585) ~ 0.6 A and EW([OIII]5007) ~ 0.5 A. In 
the case of H/3 we use Ha to identify emission as we do not expect 
to detect the former without the latter, such that A/N Ha = 4 cor- 
responds to EW(H/3) ~ 0.25 A. Thus these represent the emission 
line detection limits for our data. This method returns 1 17 galaxies 
with significant emission, 63% of the initial 186 flagged using the 
SDSS data. 

We then use GANDALF to correct, by subtracting off the best 
gaussian fit, any line with A/N ^ 4. For galaxies with detectable 
Ha emission but no significant H/3 emission we choose to force a 
H/3 correction. We have checked and found that all Ha/H/3 flux 
ratios are compatible, within errors, with Case B recombination 
dOsterbrock & Ferlano1l2006T) . During the testing of this technique 
we found an issue whereby the emission lines of certain galaxies 
were poorly fit, often, but not always, because the line in ques- 
tion was exceptionally strong, and as such left large residuals in the 
cleaned spectra. To quantify the quality of each fit we compute a 
X 2 . between the best fit and observed spectrum in the immediate 
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Figur e 2. Emission-line analysis diagrams: The left panel displays the line ratio diagnostic space of|5aldwin et al.| j 1 98 lh with the demarcation line of 
Kauf fmann etafl (2003) used to classify the 110 galaxies in our emission-line sample. Red circles are AGN/LINER while blue circles are star-forming 
galaxies. The middle and right panels show the fraction of emission line galaxies which are AGN/LINER (red lines) or star-forming (blue lines) as a 
function of luminosity and cluster-centric radius. The black dashed lines indicate the fraction of galaxies for which a definite classification could not be 
made (referred to as ambiguous emission driver, or AED, galaxies in the text). 



region around each line and visually inspect the corrected spectra 
to find the maximum acceptable In practise the most important 
line to adequately clean is H/3 as it is our primary age indicator and, 
together with [OIII]5007, one of the strongest emission lines in the 
wavelength range covered by the Lick indices. With this in mind 
we determined xlfa ^ 1.1 to provide suitable cleaning. 

To summarise, our final sample consists of 246 Coma mem- 
bers with S/N > 20 A -1 , median S/N = 35 A" 1 , and no emission 
down to the limits defined above, and 1 10 with S/N 20 A -1 , me- 
dian S/N = 34 A -1 and Xh/3 ^11 with detectable emission in any 
of H/3, [OIII]5007, Ha or [NITJ6585. We note that both spectral 
classifications, passive or emission line, are based on the spectrum 
obtained from the region of each galaxy sampled by the 3" (1.4 
kpc) SDSS fibre and may not apply more globally to the galaxy as 
a whole. 



3.2 Emission line Analysis 

Having determined which galaxies in our sample show emission 
within the spatial coverage of the SDSS fibre, we next seek to dis- 
cern what is driving said emission. To this end, in Fig. [2] we plot 
our G ANDALF derived line ratios on a version of the Baldw hTet al.l 
dl98lh or BPT diagnostic dia gram. Next we use the empirical rela- 
tion of lKauffmann et alj 42003), derived using SDSS spectroscopy, 
to delineate between star-formation and AGN/LINERs as the pri- 
mary emission driver. As not all of our galaxies have significant 
detected emission in all of the four indicators, we only attempt to 
classify those with A/N > 3 in H/3, [OIII]5007, Ha and [NII]6585 
and galaxies with A/N ^ 3 in Ha and [NII]6585 for which it is pos- 
sible to comment using [NII]/Ha alone. Applying these constraints 
we assign 51 galaxies (46% of the emission line sample defined 
previously) to be AGN/LINERs, 32 (29%) to have ongoing star- 
formation and 27 (25%) to be ambiguous based on this data (here- 
after referred to as galaxies with an ambiguous emission driver or 
AED galaxies). In Fig.|3]we show the distribution of our sample on 
the sky delineated by spectral classification. 

The left panel in Fig. [2] is qualitatively very similar to 
the equivalent di agram presented in figure 7, top left panel, of 
ISmith et al.l j2007l) who studied a large sample of galaxies in the 
Shapley supercluster, and, as with this work, imposed no selection 
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Figure 3. Distribution on the sky of the 356 galaxies in our final sam- 
ple (see text for details). Galaxies are split into passive (red filled circles), 
AGN/LINER (cyan), star forming (blue crosses) and those with an ambigu- 
ous emission driver (black points) based on the approach presented in Sec- 
tion |3.2| The large black symbols highlight the positions of NGC 4874, 
NGC 4889 and NGC 4839. 

criteria based on colour or morphology. Aside from a small number 
of star-forming galaxies toward the central left region of the BPT 
diagram in Fig. [2] the majority of cluster members with emission 
define a continuous sequence from low ratios to high ratios in bot h 
[NII]/Ha and [OIII]/H/3. As commented by ISmith et al.l fc007h . 
this makes the assignment of a star-forming or AGN/LINER clas- 
sification very sensitive to the demarcation relation employed. 

Also in Fig. [2] we plot the fraction of emission line galaxies 
which are determined to be star-forming or have AGN/LINERs rel- 
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ative to the entire emission sample, as a function of r-band absolute 
magnitude and cluster-centric radius with ~ 20 galaxies per bin. 
For reference these panels also display the fraction of AED galax- 
ies. A well defined trend is seen with luminosity in the sense that 
the fraction of AGN/LINERs increases with increasing brightness 
whereas t he fraction of star - forme rs de creases, a result in agree- 
ment with lKauffmannetal.1 d2003l) and ISmith etal] d2007l) . Inter- 
estingly no such correlation is observed as a function of radius with 
the fractional distributions remaining roughly constant, except for 
the inner-most bin where the star-forming fraction appears to di- 
minish. However, this drop coincides with a rapid increase in those 
galaxies whose emission driver is undetermined making firm con- 
clusions difficult. 



3.3 Velocity dispersions 

The SDSS spectroscopic pipeline measures velocity dispersions 
within the survey's 3" diameter fibre for the majority of spec- 
tra which are typical of early-type galaxies, and have redshifts 
z < 0.4, using a direct fitting methotff]. However, 15 (6%) of 
our non-emission sample and 36 (33%) of our galaxies found to 
have emission do not have measured velocity dispersions in DR7. 
As such we opt to use the penalised pixel fitting code (PPXF, 
Cappellari & Emsellem 2004) to obtain the velocity dispersions for 

our entire sample in a homogeneous fashion. 

Details of the technique used here can be found in lPrice et al.l 

(2009). Briefly the routine uses the same 85 stellar templates used 
by GANDALF, modulated by multiplicative and additive Legendre 
polynomials of orders 6 and 2 respectively and broadened by a 
parametric line-of-sight velocity distribution, in this case a gaus- 
sian, to construct a model of each galaxy's spectrum. Prior to fit- 
ting, the templates are matched to the variable resolution of the 
SDSS spectra. Fortunately the SDSS pipeline itself records spec- 
tral resolution as a function of wavelength for each fibre and so we 
use this information to smooth the templates with a variable width 
gaussian, typically with a lnst ~ 58-70 km s' 1 in the 4000-6000 
A fitting range used here, to match our data. Taking advantage of 
fitting in the pixel space we mask bad pixels and the NaD line at 
5892 A which can often be affected by interstellar absorption. 

To assess the uncertainty on our velocity dispersion measure- 
ments ideally one would like to undertake a bootstrap procedure, 
using the error spectra to generate new pseudo-observations and 
then passing them through PPXF, but as our sample size is large 
this approach would be time consuming. As a compromise we cre- 
ate a subsample of galaxies, 10 per S/N bin with a bin size of 10 
between S/N = 10-50, randomly selected from the full sample. For 
each selected galaxy we generate 50 realisations, run them through 
PPXF and take the standard deviation of the resulting distribution 
as the error on that galaxy's velocity dispersion. Finally we inter- 
polate the <7-S/N-<j error distribution formed by our subsample into 
a surface using a bivariate spline. It is then a simple matter of in- 
puting the S/N and a of each galaxy from the full sample to obtain 
an estimated a error . 

Comparing our velocity dispersion measurements to those 
from the SDSS pipeline for galaxies with both we find a mean offset 
<SDSS-This work> = -6 km s _1 , likely attributed to our different 
fitting range and masking procedures. Once this systematic shift is 
accounted for we find an RMS scatter of 5 km s _1 with an intrin- 
sic component of 0. Comparing our measurements with those from 



iMoore et all J2002h we find a mean offset of 2.8 km s~\ an RMS 
scatter accounting for the offset of 8 km s _1 with an intrinsic com- 
ponent of 4.7 km s" 1 . Thus, overall, we find good compatibility 
bet ween our measurem ents and those made by the SDSS pipeline 
and lMoore etafl J2002I) . 

The final step necessary before using the measured velocity 
dispersions in our analysis is to aperture correct them to a com- 
mon sampling in terms o f each galaxy's half-lig ht radius. Here we 
follow the prescription of |j0rgensen et alJ dl995). 



lo« 



&SDSS 
a R c /2 



-0.04 log 



1.5" 

Re/2 



where R e is the psf-corrected half-light or effective radius, ob- 
tained by fitting a single Sersic model using GALFIT jPeng et all 
|2002|) to the SDSS r-band image of each galaxy, usdss and <Jr s /2 
are the galaxies' velocity dispersions through the 1.5" radius SDSS 
fibre and the desired R e /2 aperture respectively. Corrections are in 
the range -0.067 to 0.017 dex with a mean of a cor = -0.004 dex or 
~ 1 %. 



3.4 Absorption line indices 

We use the LICK Ew code which is provid ed as part of the Ez- 
Ages packag^fl l lGraves & Schiavonl 12008) to measure the Lick 
in dices of our sam ple. The full s et of indices defined in Tab le 1 
of lWorthevI d 19941) and Table 1 of IWorthev & Ottavianil J 19971) are 
measured where possible, although a small fraction of the galax- 
ies have regions of bad pixels which prevent access to all indices. 
The ro utine uses the error spectra and the equations of Cardi el et al.l 
d 19981) to compute error estimates for the indices. 

Having obtained the raw indices for our sample it is then 
necessary to correct them to the resolution of the Schiavon mod- 
els, which are themselves based on stellar spectra smoothed to the 
wavelength-dependent Lick/IDS resolution. Effectively one wants 
to know the value of each observed index for a given galaxy at the 
model's resolution but without any doppler broadening from that 
galaxy's stellar velocity dispersion. One route is to employ mul- 
tiplicative and additive corrections derived from artificially broad- 
ened stars observed with the same instrument as the galaxy spec- 
tra. More recent works have moved toward using stellar population 
model templates which permit a better match in terms of spectral 
type to observed galaxy spectra, while still maintaining the zero 

velocity dispersion requirement. 

Here we opt to follow a similar method to that of lKelson et al.l 

(2006) and also take advantage of the functionality of PPXF in de- 
riving our velocity broadening corrections. First, PPXF is used to 
output the best fit combination of our stellar templates at the reso- 
lution of our data for each galaxy with and without smoothing by 
the velocity dispersion of that galaxy. Next LlCK_Ew is executed 
three times; first on the observed spectrum (I {, s ), then on the best 
fit spectrum with velocity dispersion broadening (It) and finally 
on the broadening free spectrum smoothed to t he Lick re s olutio n 
(hick) for each index as defined by Table 1 of ISchiavonl d2007l) . 
The corrected indices for each galaxy (I CO r) can then be obtained 
via: 



lobs +{h 



(1) 



3 http://www.sdss.org/dr7/algorithms/veldisp.html 



4 http://www.ucolick.org/~graves/EZ_Ages.html 
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Figure 4. Derived velocity broadening corrections for the seven indices used in our stellar population 
emission line galaxies respectively. Here the y-axis refers to the bracketed term in equation {T}. 




Figure 5. Lick indices used in our stellar population analysis, plotted against velocity dispersion for our Coma galaxies. Red filled circles are our passive galaxy 
subsample, with the four blue passive galaxies highlighted by black circles. Our emission-line sample is denoted based on emission driver as determined in 
Fig. [2] with cyan filled squares representing AGN/LINERs, blue crosses for star-forming systems and black points where the source is ambiguous. Red and 
blue error bars show typical (median) index errors for passive and emission galaxies respectively. The slope of the best fit trend line to the passive galaxy data 
is displayed in each panel where the first value includes the blue subset and the second value does not. The latter fit is also plotted as a red line in each panel. 
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In Fig. [4] we present the index corrections, that is the value of the 
bracketed term on the right hand side of equation (QJ, for the indices 
which were used in our stellar population analysis (see Section[4} 
as a function of galaxy velocity dispersion for both emission and 
non-emission line samples. Uncertainties on these corrections are 
obtained using a similar approach to that used previously for the 
velocity dispersion measurements but in practice are ^ 10% of the 
typical index errors and so are not propagated further. 

From Fig.|4]it is evident the corrections for all indices behave 
broadly as expected, that is negative corrections when a galaxy's 
velocity dispersion is sub-Lick resolution for a given index moving 
to positive corrections at super-Lick resolution. Just by inspection 
there appears to be a significantly larger scatter in the corrections 
of the emission line galaxies relative to their non-emission coun- 
terparts. This is caused by the corrections also being a function of 
index strength and, therefore, both samples in fact have comparable 
scatter. However, an increased scatter is noted at a ^ 150 km s _1 
which can likely be attributed to a greater variation in fitted spectral 
type for galaxies with lower velocity dispersion. 

As the Schiavon models are themselves based on a flux cali- 
brated stellar library, with no transformation to the Lick/IDS sys- 
tem, we also make no attempt to calibrate our data on to the Lick 
system. 

In Appendix |A] we present a comparison of our index mea- 
surements with two other literature studies and show that our data 
is generally in good agreement with both datasets. In Appendix iDl 
we tabulate our index and velocity dispersion data. 



3.5 Absorption line Analysis 

In this section we shall study the stellar populations of our sam- 
ple galaxies in a model independent way, using their absorption 
line indices. We will look for correlations between these measure- 
ments and galaxy velocity dispersion and local environment, pa- 
rameterised by cluster-centric radial distance. The same analysis 
techniques are also applied in Sections 14.41 and l4~5l to the model 
dependent stellar populations parameters. 

In Fig. [5] we plot the indices used in our analysis (where 
<Fe>=(Fe5270+Fe5335)/2) against velocity dispersion and in 
Fig. [6] against cluster-centric radius. In the former figure the leg- 
end of each panel displays the results of a linear least-squares fit 
to the passive galaxy data weighted by the index uncertainties and 
iteratively taking into account intrinsic scatter. Slope uncertainties 
are computed by creating 500 realisations of the index data using 
their errors and taking the standard deviation of the resulting dis- 
tribution of slopes. The first reported slope is computed including 
the blue subset identified earlier while the latter slope is derived 
having masked them out, resulting in a fit which is more com- 
parable to previous works based on colour-selected samples. For 
the index-radius panels the index data are binned for clarity due 
to substantial intrinsic scatter and the emission-line sample is split 
by emission driver. The passive sample has bins derived excluding 
the four blue quiescent outliers and contains ~ 27 galaxies per bin 
while the emission line bins contain ~ 13 and ~ 7 galaxies per 
bin for the AGN/LINER and AED samples respectively. The fits to 
the star forming sample are not plotted in Fig. [6] for clarity. This 
is because they often have trends with strongly offset intercepts or 
slopes in the opposite sense to the other subsets and thus act to dis- 
tract the eye from the passive sample which is more so the focus 
of this work. The results of the index-radius fits, including the star 
formers, are displayed in Table[2] We stress that in all cases it is the 



unbinned data that are fit and the binning is for display purposes 
only. 

The relations for the non-emission line galaxies presented in 
Fig.[5]are broadly c onsistent with previou s studies on red sequence 
galaxies in Coma jMafkovic et alj[2009h . and other low redshift 
clusters dSmith et alj|2007l> . Both of these works have spectra de- 
rived from comparable physical scales to this work, the former ob- 
serving with 3.1" diameter fibres and the latter with 2" diameter 
fibres at z = 0.048, giving a spatial coverage of ~ 1.8 kpc. Strong 
positive correlations are particularly evident for CN2, C4668 and 
Mgb while a significant anticorrelation is present for H/3. There is 
evidence, at least by eye, for increased scatter about the fitted rela- 
tions with decreasing a for all indices except perhaps CN2. 

The emission-line galaxies generally have higher H/3 and 
weaker metal lines at low a with respect to the quiescent galax- 
ies while at high a both subsamples have similar index strengths. 
As displayed in Fig. [5] this trend is clearly driven by a transition 
from AGN/LINER dominated systems at high a, with stellar popu- 
lations at least comparable to quiescent bright early-type galaxies, 
to star-forming galaxies at low a. Over the range in velocity disper- 
sion sampled here the transition is fairly smooth, with a crossover 
at ~ 100 km s _1 and a marked overlap between the two emission 
drivers, a feature that acts to reiterate Fig.[2]middle panel. 

In Fig. [6] significant radial trends are seen for the passive 
galaxies in CN2 (63a), C4668 (2.8a), H/3 (4.1a) and Mgb (7.7a) 
which typically increase in significance when masking the blue 
subset. For H/3 and C4668 the correlations are positive with larger 
index values at greater cluster-centric radii whereas CN2 and Mgb 
are observed to decrease with radius. Strong correlations are also 
identified in these indices for the AGN/LINER galaxies in the same 
sense as the passive sample, with the exception of C4668. The in- 
dices of star-forming galaxies behave somewhat differently with 
H/3 decreasing and Mgb increasing with radius (see Table |2), al- 
though the robustness of these trends given the small size of the 
subsample is questionable. 

The trends presented in Fig. [5] and [6] may be affected by cor- 
relations between radius and velocity dispersion and thus, for in- 
stance, if the outer parts of the cluster are preferentially inhabited 
by lower a galaxies the derived index-radius fits will be biased. 
A further factor we seek to address stems from the fact that the 
survey's 3" fibres sample a different fraction of each galaxy based 
on their physical size. As such, if we assume some typical inter- 
nal radial gradient in their stellar populations, and therefore index 
strengths, larger galaxies will generally have index strengths bi- 
ased higher, in the case of negative gradients, or lower, in the case 
of positive internal trends, relative to smaller galaxies. Indeed, on 
average, negative gradients have been detected in the metallic in- 
dices while H/3 gradients may be positive or ne gative for individ- 
ual galaxies, but are typically nea rly flat (e.g. iDavies et al Jl 19931 : 
lGonzaled[i993l ; lRawle et all2008h . 

To factor out the issues discussed above and derive index re- 
lations at fixed velocity disp ersion and radius that are independent 
of aperture effects we follow lSmifh et al. (2006, hereafter S06) and 
S08 and fit our data with a multi-dimensional plane that incorpo- 
rates all three terms, 



I = aR cc + b log a + c log h d (2) 

tie 

where 1.5" and i? e are again the fibre radius and effective radius of 
each galaxy respectively. Here R cc is the projected cluster-centric 
radius. The coefficients of the first two terms give the desired in- 
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Figure 6. Lick indices against cluster-centric radius. Data are binned with roughly the same number of galaxies per bin. Red (~ 27 galaxies per bin) and 
cyan (~ 13 galaxies per bin) filled circles denote passive and AGN/LINER galaxies respectively. The black points (~ 7 galaxies per bin) are galaxies whose 
emission driver could not be determined. Trend lines are fit to the unbinned data. The star forming subsample is not plotted for clarity (see text). 




Figure 7. Index-radius diagrams plotted in the form of a residual index 51 where 81 = I — (felogcr + clog(1.5"/-Re) + d). The gradient of the best fit 
lines corresponds to the coefficient a in a planar fit of the form I = aR cc + b log cr + c log(1.5"/-Re) + d . Thus these diagrams compare the data at fixed 
velocity dispersion. Symbols and binning are as in Fig.|6]and once more the star forming subsample is not plotted for clarity. 



dex trends. Note, while S08 use luminosity in their plane fits, we in a comparable manner to Fig. [6] we plot the data in the form of a 
almost always find stronger index correlation s with a and therefore residual index SI, 
opt to employ the latter in our fits (see also iBernardi et al.ll20o3 : 
Smit h et al . 2009b). In order to display the results of the planar fits 

1 5" 

81 = I - (bloga + clog — — + d) 
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Table 2. Cluster-centric index gradients in units of A, or mags for CN2, per degree obtained from the fits in Fig.[6]and Fig. [7] For our passive galaxies 
two slopes are given, the first including the blue subset outlined earlier and the second masking them. The 2nd column indicates the number of 
galaxies in each respective passive galaxy fit out of a potential of 246. The number of galaxies in each emission driver delineated fit is shown in the 
header row and remains constant from one index to another. We express the significance of each fit, in units of standard error, in parentheses. 



Index 


N 


Passive 


AGN/LINER (51) 


SF (32) 


AED (27) 








I — R C C 


(Rg.|6) 






CN2 


244/240 


-Q.009±Q.001(6.3) 


-0.012±0.001(8.1) 


-0.022±0.003(8.8) 


0.009±0.005(2.0) 


-0.0 15 ±0.003(5.1) 


Ca4227 


246/242 


-0.025±0.019(1.3) 


-0.029±0.019(1.5) 


-0.061±O.036(1.7) 


0.058±0.078(0.7) 


-0.100±0.050(2.0) 


C4668 


245/242 


0.183±0.065(2.8) 


0.161±0.061(2.6) 


-0.432±0. 109(4.0) 


0.766±0.250(3.1) 


-0.263±0.140(1.9) 


Hbeta 


246/242 


0.104±0.026(4.1) 


0.132±0.025(5.3) 


0.233±0.046(5.1) 


-0.344±0. 107(3.2) 


0.234±0.061(3.8) 


Mgb 


246/242 


-0.223±0.029(7.7) 


-0.251±0.028(8.9) 


-0.289±0.044(6.5) 


0.465±0. 114(4.1) 


-0.299±0.067(4.4) 


<Fe> 


245/241 


0.023±0.025(0.9) 


0.018±0.023(0.8) 


-0.064±0.039(1.6) 


0.212±0.096(2.2) 


-0.064±0.058(1.1) 








<5I - R cc 


(Fig.EJ 






CN2 


244/240 


-0.001±0.001(0.4) 


-0.002±0.001(1.6) 


-0.008±0.002(3.4) 


0.008±().005(1.5) 


-0.006±0.004(1.6) 


Ca4227 


246/242 


-0.022±0.02 1(1.0) 


-0.025±0.020(1.3) 


-0.052±0.034(1.5) 


0.016±0.084(0.2) 


-0.095±0.059(1.6) 


C4668 


245/242 


0.345±0.060(5.8) 


0.334±0.058(5.7) 


-0.148±0.107(1.4) 


0.692±0.253(2.7) 


-0.239±0.188(1.3) 


Hbeta 


246/242 


0.044±0.026(1.7) 


0.077±0.025(3.1) 


0.118±0.045(2.6) 


-0.327±0. 114(2.9) 


0.237±0.083(2.9) 


Mgb 


246/242 


-0.099±0.027(3.7) 


-0.114±0.027(4.3) 


-0.108±0.046(2.4) 


0.467±0. 116(4.0) 


-0.216±0.082(2.6) 


<Fe> 


246/242 


0.041±0.029(1.4) 


0.035±0.023(1.5) 


-0.042±0.056(0.7) 


0.203±0.129(1.6) 


-0.048±0.093(0.5) 



Table 3. Comparison of the index gradients derived by this work (from the planar fitting technique, second section Table [2} with those reported in 
the literature for Coma by C02 and S08 and for a large number of nearby clusters by S06. For the C02 results their Mg2 gradient is converted to an 
equivalent gradient in Mgb. The results of S06 are converted to gradients in degrees assuming R200 = 1-59° for Coma. "Core-SW" denotes our sample 
restricted to these areas (see text). 



Index 


This work 


Core-SW 


C02 


S08 


S06 


CN2 


-0.002±0.002 


-0.006±0.003 




-0.025 ±0.007 




Ca4227 


-0.025±0.020 


-0.060±0.043 




-0.136±0.072 




C4668 


+0.334±0.058 


+0.622±0.150 




-0.207±0.270 


-0.102±0.046 


H/3 


+0.077±0.025 


+0.119±0.056 


+0.327±0.116 


+0.539±0.104 


+0.046±0.013 


Mgb 


-0.114±0.027 


-0.070±0.058 


~0.491±0.125 


-0.490±0.145 


-0.128±0.016 


<Fe> 


+0.035±0.023 


+0.095 ±0.052 


-0.082±0.084 


+0.031 ±0.083 


+0.()18±0.015 



We choose not to plot the equivalent 81 for b log a since the slopes 
typically change by less than one standard error relative to those 
shown in Fig. [5] However, the effect on the radius relations is gen- 
erally more substantial and so in Fig. [7] we present SI-K CC for our 
sample galaxies with fit results again tabulated in Table [2] The ro- 
bustness of our fitting procedure is examined in Appendix Icl 

For the passive galaxies significant correlations are still found 
for C4668 (5.8a), Mgb (3.7a) and H/3 (3.1cr), when the blue subset 
is masked, while radial trends in CN2 become insignificant. On the 
other hand, our sample of AGN/LINERs now generally have much 
weaker trends relative to their linear fits, with significant correla- 
tions still found for CN2 (3.4a), n(3 (2.6a) and Mgb (2.4a). This 
sharp reduction in the computed gradients stems from the fact that 
velocity dispersion correlates mildly with radius for this sample in 
the sense that lower dispersions are found at larger radii. Finally, 
the index trends for star-forming galaxies remain comparatively un- 
changed. 

A valid question here is to what extent do the outermost points 
drive the significant radial trends recovered by both the linear and 
planar fits for our passive sample. To test this we have repeated 
the fitting procedure with the outermost points, those in the final 
bin with R cc > 1.55°, and the blue subset masked. There is some 
change to the fitted slopes reported in Table [2] However, the three 
indices with significant radial coefficients in the planar fits remain 
C4668 (6a), H/3 (2.5a) and Mgb (3.1cr) and with the same sign as 
given in Table [2] Indeed, the same is true for the linear fits with 
C4668 (4.4a), H/3 (3.3a) and Mgb (5.7a). 

In Table [3] we compare the radial gradients for our passive 
sample to those reported in C02 and S08 for Coma passive galax- 
ies and more generally in a sample of low redshift clusters from 



S06. The sign of our significantly detected slopes are broadly con- 
sistent with those found elsewhere, with the exception of C4668, 
and indeed we also find at most a very weak correlation for <Fe>. 
Interestingly C02 and S08 find much stronger gradients in H/3 and 
Mgb relative to both this work and S06. This result may either be 
driven by the luminosity range covered here compared to the pre- 
vious works, C02 analysed both dwarf and giant galaxy regimes 
while S08 only include dwarf galaxies, or the spatial coverage of 
the cluster provided by both previous studies which was limited to 
core and south-west regions, or some combination of both factors. 
To try to assess the impact of limited spatial coverage we create a 
subset of our passive sample with a highly comparable cluster foot- 
print to that of S08 and refit our planar relations. Our Core-SW sub- 
set has 1 16 galaxies, 3 of which are classified here as blue passive 
galaxies. The results of these fits, with the blue galaxies masked, 
are presented in Tableland do not display a significant increase in 
H/3 or Mgb gradient, although we do note that when including the 
blue subset in the fit the derived H/3 slope increases to 0.194±0.061 
A deg -1 . This finding seems to imply that the weaker slopes found 
here for H/3 and Mgb compared to S08 are driven by our luminosity 
range being constrained to bright Coma galaxies not by our clus- 
ter coverage. Similarly, C02 including a sample of dwarfs in their 
data likely results in their slopes being an intermediate between this 
work and S08. 



4 STELLAR POPULATION PARAMETERS 

Having measured the absorption indices for our passive and 
emission-line samples we are now able to convert these into es- 
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timates of their stellar population parameters using the Schiavon 
models and Ez-AGES fitting code. As we plan to determine these 
parameters for both of our subsamples it is important to preface this 
section with a word of warning. 

It has been known for some time that the presence of compos- 
ite stellar populations (CSP) resulting from multiple star-forming 
bursts significantly affects how the s tellar age returned b y SSP 
models should be interpreted. Recentlv lSerra &Tragerl 12007) con- 
ducted extensive simulations to fully assess the impact of trying 
to apply SSP models to CSPs, which they constructed based on 
a two burst scenario with a mass dominating old population and 
a frosting of younger stars. They found that rather than obtain- 
ing a luminosity-weighted age, the SSP-equivalent age is often 
skewed substantially lower by a relatively small mass fraction of 
younger stars, an outcome attributed t o the non-linear response of 
the Balmer lines to stellar age (see also lTrager & Somervillel2 009 , 
for a more detailed discussion). However, they report that SSP- 
equivalent [Fe/H] and [a/Fe] track quite closely the luminosity- 
weighted parameters, being driven typically by the older popula- 
tion except in the most extreme cases of a very young and mas- 
sive frosting. We therefore stress that, particularly in the case of 
the star-forming galaxies in our emission-line sample, the derived 
ages are Balmer- weighted SSP-equivalent parameters and may not 
be indicative of the global stellar populations in these galaxies. 

4.1 The Model and Parameter Fitting 

The Schiavon models are based on the |jonesl ll999) empirical stel- 
larlibrary and span a range of 1 to 17.7 Gyr and -1.3 to +0.2 [Fe/H]. 
They enable a determination of abundance ratios in terms of Mg, 
C, N and Ca by incorporating their stellar atmospheric effects. As 
the models are based on scaled-solar isochrones the evolutionary 
tracks and stellar atmospheres are not treate d in a self-consistent 
mann er in a similar vein to the widely used iThomas et alj d2003l . 
120041 hereafter TMBK) models. Here we employ a Salpeter initial 
mass function. 

We derive the SSP-e quivalent parameters for ou r galaxies us- 
ing the Ez-AGES code of iGraves & Schiavon! d2008h which, using 
the Schiavon models, seeks to return consistent stellar population 
parameters across a number of index-index diagrams. The routine 
begins by making an initial estimate of age and [Fe/H] by invert- 
ing a H/3 - <Fe> grid, specifically chosen to strongly break the 
age-metallicity degenera cy while being largely insensitive to non- 
solar abundance patterns dSchiavonl2 007). Next the code moves on 
to determining [Mg/Fe] from a H/3 - Mgb grid. To do this it first 
inverts the new grid, compares the age and [Fe/H] measured from 
this grid to the fiducial estimates and if they do not agree assumes a 
non-solar [Mg/Fe] is present and recomputes the model with a new 
[Mg/Fe]. Further iterations occur until the derived age and [Fe/H] 
converge to the fiducial values within some tolerance. The code 
then goes on to repeat this iterative approach replacing Mg by C, 
N and Ca in turn with the final abundance pattern used to compute 
a new model and derive a new fiducial age and [Fe/H], with the 
entire fitting process repeating until the fit does not improve. The 
standard index set used for the fitting is H/3 and <Fe> for the ini- 
tial age and [Fe/H] estimates, Mgb for [Mg/Fe], C4668 for [C/Fe], 
CN2 for [N/Fe] and Ca4227 for [Ca/Fe]. 

We could not obtain fits for 1 1 % of our entire sample as they 
fell off the fiducial H/3 - <Fe> grid. This may be because they 
either have SSP parameters, perhaps only marginally, outside the 
models coverage or else were carried off the grid by measurement 
errors. As a work around we fit these galaxies using Fe5270 or 
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Figure 8. Histogram showing the distribution of uncertainties on our de- 
rived SSP-equivalent ages for galaxies in our quiescent sample. 

Fe5335 only instead of their mean. Typically the results obtained 
from the three different approaches are consistent, within errors, 
for those galaxies which can be fit by all three. 

To compute uncertainties on the stellar population parameters 
we select five quiescent galaxies from our original sample of 417 
with median S/N per A ~ 15, 25, 35, 45, 55 which have index er- 
rors representative of the entire sample in their respective S/N bin. 
Next for each S/N step we generate 300 new realisations of the in- 
dex data by randomly perturbing the initial values by their errors 
and run these through Ez-AGES. At each S/N we then assign the 
typical uncertainties for age, [Fe/H] and the abundance ratios as the 
standard deviation of the distribution of parameters obtained from 
our simulations. Finally we linearly fit to these points a relation of 
the form P RTr0 r = m(S/N)~ 1 + c, where P e „-or is the standard 
error on P and P = log(age), [Fe/H], [Mg/Fe], [C/Fe], [N/Fe] and 
[Ca/Fe] in turn. Therefore using this error model and a galaxy's S/N 
per A in the wavelength range 4000-6000A we are able to estimate 
the statistical uncertainty on all its stellar populations parameters. 
This approach also acts to assess the error correlations between fit- 
ted parameters. As an example, in Fig. [8] we show the distribution 
of the standard error on log(age) for our quiescent sample. This 
figure demonstrates that our S/N cut results in at most a ~ 50% 
uncertainty on our age estimates, with a median error of ~ 40%. 
Note here we assume that the index error distributions for the qui- 
escent and emission-line sample are comparable and as such do 
not attempt to quantify any systematics generated by our emission 
cleaning technique. 

4.2 Identifying Composite Stellar Populations 

As outlined earlier, the measured stellar age reported by our analy- 
sis can be significantly biased by a frosting of hot young stars and 
in such cases is far from a luminosity-weighted parameter. Fortu- 
nately tools are available that allow the detection of particularly 
young (< 1 Gyr) star bursts, thus placing a lower limit on the ces- 
sation of star formation in a given galaxy. This is particularly appli- 
cable to our passive sample since we wish to test that if a relation 
exists between age and velocity dispersion it is not driven by, for in- 
stance, the ages of lower a quiescent galaxies being biased younger 
by recent star formation. 



12 J. Price et al. 




logo- 



Figure 9. The difference between SSP-equivalent ages derived using H7 
and H/3 indices plotted against galaxy velocity dispersion for our passive 
sample. The solid line represents a simple linear fit to the data, the slope 
of which is only mildly significant (~ 2tr). The panel indicates that it is in 
fact high a galaxies which typically have a weak discrepancy between their 
Balmer line ages rather than galaxies with low velocity dispersions. 
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Figure 10. Rose Call index against H/3 for passive sample galaxies with er- 
rors in Call ^ 0.15. The parameter coverage of simulated CSPs with frost- 
ing mass ^ 1% and age ^ 0.7 Gyr and frosting mass ^ 3% and age ^1.0 
Gyr are denoted by blue and green shaded polygons respectively. The re- 
gion spanned by genuine SSPs with 2 Gyr < age ^ 14 Gyr is highlighted 
by a red shaded polygon. Despite sizable scatter in Call, the majority of our 
passive sample show little evidence for episodes of star formation within 
the past ~ 1 Gyr. 

One such method involves the comparison of SSP-equivalent 
ages derived using different Balmer indices, relying on the fact 
that the higher order Balmer lines are more strongly affected by 
the presence of a young, blue subpopulation. Therefore, using US 
or H7 rather than H/3 as our age sensitive index should result 
in younger ages being obtained if a sufficiently young frosting is 
present. Here we employ H7 derived ages which are measured by 
Ez-AGES following the conclusion of abundance pattern fitting. In 
Fig. [9] we plot the difference between the H7 and H/3 based age 



Table 4. Measured stellar population parameters for all subsamples, includ- 
ing those for early-type galaxies in the passive sample based on morpholo- 
gies taken from NED. Columns 2-5 give the parameter means (p), total 
scatter about the mean (crtot), scatter due to errors (cr errs ) and intrinsic 
scatter in the distribution (<Tj n t) respectively. The number of galaxies in 



each subset is given in parentheses. 


Parameter 


M 


°tot 






Passive (222) 


log(age) 


0.85±0.02 


0.25 


0.15 


0.20 


[Fe/H] 


-0.08±0.01 


0.17 


0.13 


0.11 


[Mg/Fe] 


0.15±0.01 


0.11 


0.11 


0.00 


[C/Fe] 


0.15±0.01 


0.15 


0.10 


0.11 


[N/Fe] 


0.04±0.01 


0.19 


0.12 


0.14 


[Ca/Fe] 


0.00±0.01 


0.12 


0.14 


0.00 


Passive Early-type (160) 


log(age) 


0.88±0.02 


0.25 


0.15 


0.20 


[Fe/H] 


-0.08±0.01 


0.17 


0.12 


0.11 


[Mg/Fe] 


0.17±0.01 


0.10 


0.11 


0.00 


[C/Fe] 


0.15±0.01 


0.16 


0.10 


0.12 


[N/Fe] 


0.05±0.01 


0.19 


0.12 


0.15 


[Ca/Fe] 


0.00±0.01 


0.12 


0.14 


0.00 


AGN/LINER (44) 


log(age) 


0.70±0.05 


0.31 


0.15 


0.27 


[Fe/H] 


-0.03±0.03 


0.18 


0.12 


0.13 


[Mg/Fe] 


0.16±0.02 


0.10 


0.10 


0.02 


[C/Fe] 


0.14±0.03 


0.17 


0.10 


0.15 


[N/Fe] 


0.06±0.02 


0.16 


0.11 


0.12 


[Ca/Fe] 


-0.04±0.02 


0.13 


0.13 


0.04 


SF (24) 


log(age) 


0.26±0.03 


0.14 


0.17 


0.00 


[Fe/H] 


-0.52±0.08 


0.41 


0.16 


0.37 


[Mg/Fe] 


0.13±0.04 


0.19 


0.15 


0.12 


[C/Fe] 


-0.05±0.04 


0.18 


0.13 


0.12 


[N/Fe] 


-0.17±0.06 


0.31 


0.16 


0.26 


[Ca/Fe] 


-0.08±0.05 


0.24 


0.20 


0.14 


AED (26) 


log(age) 


0.78±0.04 


0.21 


0.14 


0.16 


[Fe/H] 


-O.03±0.03 


0.14 


0.11 


0.08 


[Mg/Fe] 


0.19±0.02 


0.09 


0.10 


0.00 


[C/Fe] 


0.21 ±0.03 


0.15 


0.09 


0.12 


[N/Fe] 


0.08±0.04 


0.19 


0.11 


0.15 


[Ca/Fe] 


0.00±0.02 


0.09 


0.12 


0.00 



estimates for the 205 passive galaxies with both, as a function of 
their velocity dispersion. The diagram demonstrates that the ages 
derived for low a galaxies are consistent, despite sizable scatter, 
and that it is in fact the high velocity dispersion galaxies which 
show a mild discrepancy, having H7 ages at most some ~ 25% 
younger than those from H/3. Furthermore, the distribution of age 
differences about zero is dominated by a scatter due to measure- 
ment errors of 0.18, with a small in trinsic scatter of 0.07. 

Following ISmith et al.l d2009bl) , a second technique to search 
for the pre sence o f rece nt star bursts in our passive sample makes 
use of the iRosd 1 19841) Call index. This index is computed as 
the ratio of the residual flux at the core of the Ca K and H lines 
and is approximately constant for stars of later type than F2, 
but changes strongly for hotter stars. Thus, when combined with 
Balmer line indices, this index can be used to distinguish near 
SSP-like galaxies from those with a small frostin g of younger stars 
plus a mass-dominant underlying old population dLeonardi & Rose! 
1996). Indeed, such a test also places limits on the contribution 
fro m any other A-type stars, s uch as hot horizontal branch stars 
(e.g. lMaraston & Thomasll2000l) . 

Taking this approach, in Fig.[l0]we plot Call against H/3 for 
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the 149 (67%) of our passive sample fit by Ez-Ages (see next sec- 
tion) with uncertainties in their Call index = SC 0.15. Next we simu- 
late a set of CSPs using the SSP models of IVazdekis et alj J2010h 
with a base population of age > 5 Gyr and [Fe/H] > -0.68 and frost- 
ings of age ^ 1 Gyr and [Fe/H] > -0.68. The frostings have mass 
fractions of 1%, 3% and 5%. The measured Call and H/3 indices 
for all combinations of frosting mass ^ 1% and age ^ 0.7 Gyr and 
frosting mass ^ 3% and age ^ 1.0 Gyr are plotted on Fig.QJJ] For 
clarity the regions in the parameter space spanned by these models 
are represented by blue and green shaded polygons respectively. Fi- 
nally we also plot the coverage predicted for genuine SSPs with 2 
Gyr < age ^ 14 Gyr and -0.68 < [Fe/H] ^ 0.2, roughly compara- 
ble to the parameters observed for our Coma passive sample (see 
Fig.HU, in red. Of course here we note the mixing of two model 
sets, Vazdekis and Schiavon, although we expect the effects to be 
limited. 

Fig.llOlshows that, despite sizable scatter in Call, the majority 
of the 149 galaxies tested do not show evidence for small amounts 
of recent star formation. Indeed, by resampling the data given the 
measurement uncertainties we confirm that 88±2% have Call ^ 
1.1 or H/3 ^ 2 and thus are outside the parameter space covered 
by the explored frosting options. Furthermore, these galaxies have 
<CaII> = 1.150 ± 0.01 which may be compared to <CaII> = 
1.140 for the plotted SSP models. That said, ~ 15% of the sample 
have 2.0 ^ H/3 ^ 2.5 and could reach the blue and green polygons 
within their la errors. However, it is important to note that, while 
frosting may be present to some extent in these systems, many of 
them are also consistent with SSPs within their uncertainties. For 
reference, the most extreme galaxy on this plot, H/3 ~ 4.5 and Call 
~ 0.7, is the k+a galaxy GMP 3892. 

To summarise, we have found little evidence that the majority 
of the passive sample are having their SSP-equivalant ages biased 
low by small mass fractions of young stars from recent star for- 
mation episodes. We cannot, however, exclude this scenario for a 
small fraction of the sample. 

4.3 Results 

In Fig.QT]we plot the index data for both our passive and emission- 
line samples, again separated by emission driver, on an example H/3 
- <Fe> model grid. This diagram contains the 222 (90%) passive 
and 94 (85%) emission-line galaxies for which Ez-AGES yielded 
results. The vast majority of the failed fits continue to be caused by 
galaxies falling off the fiducial Balmer grid, regardless of which of 
the three Iron indices (Fe5270, Fe5335 or the mean of the two) were 
used. Also, we were unable to measure one or more of the required 
indices for a full abundance pattern fit for four of our passive sam- 
ple due to bad pixels in the respective index bandpasses. We note 
that the model grid presented in Fig. QT] is for a solar-abundance 
model only but, as this index pair is specifically selected for its 
weak sensitivity to non-solar elemental abundances, this particular 
grid is relatively stable in index space. 

Simple inspection of Fig. QTJ indicates that quiescent and 
AGN/LINER dominated galaxies have a comparable spread in age 
and [Fe/H] which, from visual assessment alone, typically spans 
from age> 3 Gyr and [Fe/H] > -0.7 to the edge of the model 
grid. On the other hand the star-forming subsample are generally 
younger, age^ 3 Gyr, with a spread in metallicity that covers the en- 
tire grid, -1.3 < [Fe/H] < +0.2. In Table|4]we compare the subsam- 
ples in a more quantitative manner, presenting the mean, rms scat- 
ter (atot), scatter due to errors (a errs ) and intrinsic scatter (oi nt ) 
of the six stellar population parameters measured here. The latter is 
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Figure 11. H/3 and <Fe> indices for the 316 (89%) of sample galaxies 
which were fit by Ez-AOES overplotted on a Schiavon model grid with 
a solar abundance pattern. Red filled circles, Cyan filled squares and blue 
crosses are passive, AGN/LINER and star forming galaxies respectively. 
Black points represent galaxies whose emission driver could not be deter- 
mined. Red and blue error bars show median index errors (oup = 0-21 
and cr<Fe> = 0.19 in both cases) for passive and emission line galaxies 
respectively. Due to the relative insensitivity of this particular index pair to 
abundance pattern variations the grid is fairly stable in index space. 



computed as af nt = crf ot - a1 rrs . A number of noteworthy points 
emerge from this table. 

Firstly, our passive galaxy sample shows a large age spread 
of ±0.19 dex (3.1 Gyr) even after factoring out scatter due to sta- 
tistical error. While not directly comparable to the recent work of 
T08 who selected their sample based on morphology, this large in- 
trinsic scatter in the ages of Coma passive galaxies, the vast ma- 
jority of which occupy the cluster's red sequence, does not sup- 
port their finding of age uniformity across the range in luminosity 
and a studied here. They also conclude that differing sample se- 
lection techniques in terms of morphology (early-types) vs colour 
(red sequence) are unlikely to be the main driver behind this dis- 
crepancy. To test this we obtained morphologies from NED for as 
many galaxies in our passive sample as possible and find that ~ 
72% of the 222 galaxies are of types E-S0. The mean stellar popu- 
lation parameters for the passive early-types are presented in Table 
[4] Since the two samples are not independent, we use a resampling 
method to assess the significance of the differences between them. 
Only the mean log(age) and [Mg/Fe] are found to be different at the 
~ 3a level. This test also confirms that the observed intrinsic age 
scatter is not driven by our sample including emission-free spiral 
galaxies. 

Secondly, the AGN/LINER systems have a younger mean age 
than the passive galaxies yet similar mean [ Fe/H] and elementa l 
abundanc es, in agreement with the findings of lGraves et al. N2007b . 
Of course iGraves et at] J2007t) apply a colour constraint when se- 
lecting their sample and so, as our analysis is based typically on 
red sequence quiescent galaxies and AGN/LINERs which display 
a range of colours, a direct comparison may be inappropriate. To 
test this we apply the same colour cut used to identify the blue 
passive subset and find that, for the 40 galaxies that meet this cri- 
teria, their mean stellar population parameters remain effectively 
constant with fJ,i os ( age ) changing by less than one standard error 
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(Miog(age) = 0.72 ± 0.05 dex). Thus AGN/LINER galaxies that are 
on the cluster's red sequence are ~ 15-20 % younger than our pas- 
sive sample. One further caveat here, however, is that the two sub- 
samples may have different velocity dispersion distributions which 
perhaps introduces a further bias. This issue will be addressed in 
the next section. 

Thirdly, the 26 Coma galaxies with unclassified emission have 
average stellar population parameters similar to those of the pas- 
sive and AGN/LINER samples. Indeed they are found to have a 
mean age between that of the former and latter. This likely rules 
out the presence of weak ongoing star formation in these galax- 
ies and points toward low-level AGN/LI NER activity. This p redic- 
tion is further reinforced by the results of lGraves~ et al. ( 2007]) who 
find that weak AGN/LINERs display ages closer to those of pas- 
sive galaxies (see their Fig. 12, top row of panels) and that, while 
not plotted in Fig. [2] due to low line A/N, the AED galaxies all 
have [NII]/Ha and [OIII]/H^ ratios compatible with AGN/LINER 
activity. 

4.4 Trends with Velocity Dispersion 

In this and the following section we seek to identify if the stellar 
population parameters of our sample galaxies display correlations 
with velocity dispersion and cluster-centric radius. To this end we 
will conduct a similar analysis to that carried out for the absorption 
line indices, fitting for both linear and planar relations. 

In Fig. [12] we present the six measured stellar population pa- 
rameters as a function of galaxy velocity dispersion for both pas- 
sive and emission line samples. The legend of each panel again 
displays the results of a weighted linear fit to, firstly, the 222 galax- 
ies in the quiescent sample and, secondly, having masked the two 
blue passive galaxies that were fit by Ez-AGES. Once more these 
fits are obtained taking into account intrinsic scatter with the slope 
uncertainties being computed by resampling the stellar population 
parameter data using their errors. As with the index-er diagrams 
presented earlier, we find that the velocity dispersion coefficient in 
a planar fit is always well within one standard error of the slope 
derived from its respective linear fit. We therefore choose to only 
plot the linear relations but note that the discussion that follows also 
applies to trends corrected to fixed radius and taking into account 
aperture effects. 

The first panel of Fig. [T2] indicates a strong positive correla- 
tion between SSP-equivalent age and velocity dispersion for Coma 
passive galaxies. From visual inspection it is apparent that the scat- 
ter about the fitted relation increases for lower a galaxies to the 
extent that a fraction of these systems have ages comparable to 
galaxies with the highest velocity dispersions in the cluster. Indeed, 
low a galaxies are seen to span the full range of ages provided by 
the Schiavon models. Extensive scatter in the ages of Coma low a 
galaxies has also be en report ed by T08, based on index data from 
iMoore et al.l J2002h . and by ISmith et al.l J2009ij) for their dwarf 
sample. This implies that passive galaxies with low velocity dis- 
persions in Coma show a wider range of star formation histories 
than do high a galaxies. 

This result is unlikely to be driven by an increased incidence of 
star bursts with ages < 1 Gyr in low a systems as confirmed earlier. 
As before, fitting separately to the early-type subsample does not 
significantly affect the observed relation, in fact we find a slope 
that is highly consistent with that obtained for the entire passive 
sample. Furthermore, the trend is altered by less than one standard 
error when derived from a planar fit. 

One additional factor deserves further comment. Kels on et al.l 



J2006h pointed out that, at fixed <r, the stellar population parameters 
of a magnitude-limited sample will be biased in the sense that, for 
example, older galaxies are fainter and so may be excluded by the 
limit. This issue will of course be more prevalent at lower veloc- 
ity dispersions where an increasingly larger fraction of old systems 
are too faint and fall out of the sample. In the first panel of Fig. 
[12] this would then result in the fitted slop e being biased too s teep. 
However, it has been d emonstrated by both l Smith et alJ ( 2009b) and 
lAUanson et alj 1 120091) that such a selection effect does not strongly 
influence simple two-parameter fits and so it is unlikely to signifi- 
cantly alter the recovered relation. 

Having considered all the points discussed above, we must 
conclude that the data presented here support a correlation be- 
tween galaxy SSP-equivalent age and velocity dispersion, or in 
other words archaeological downsizing, for passive galaxies in the 
Coma cluster across the studied a range. 

In this panel we also fit separately to the red AGN/LINER 
galaxies identified previously, denoted by the blue dashed line in 
Fig. [T2] the result of which clearly demonstrates that at fixed ve- 
locity dispersion these systems are, on average, younger than their 
passive counterparts. 

The second panel of Fig. [12] displays a much weaker positive 
correlation, relative to the age-cr relation, between [Fe/H] and a 
for the passive galaxies. There is some evidence for an asymmet- 
ric scatter toward lower [Fe/H] at the low a end of the relation, 
although we note that the upper bound is inherently limited by the 
models coverage which is capped at [Fe/H] = 0.2. AGN/LINER 
systems appear to be distributed in a similar manner to the pas- 
sive galaxies in this parameter space while the star forming sam- 
ple shows only a limited overlap, with a sizeable fraction hav- 
ing [Fe/H] ^ -0.5. In fact a number of the latter sample have 
[Fe/H] < -1, a regime where Ez-AGES, and presumably the Schi- 
avon models, may well be unsuitable as cautioned by the authors 
dGraves & Schiavorj|2008l) . 

In the third panel of Fig. [T2] we see that higher a quiescent 
galaxies tend to have higher [Mg/Fe], and therefore shorter star 
formation timescales, than lower a galaxies. The fitted slope in thi s 
panel is highly comparable to that derived bv lSmith et al.l l2009al) , 
an interesting result given t he different velocity dispersion cover- 
age of the two samples, and lGraves et al.l (2007J) for giant galaxies 
based on stacked SDSS spectra. The majority of the AGN/LINER 
and star forming galaxies scatter across a similar range in [Mg/Fe] 
to the passive galaxies with only one sizeable outlier, likely ex- 
plained by this galaxy also having [Fe/H] = -1.18. 

The final three panels of Fig. [12] show relatively strong pos- 
itive correlations for [C/Fe]-<j and [N/Fe]-cr while [Ca/Fe] is seen 
to only weakly correlate with velocity dispersion over the range 
studied here. Once more, the emission line galaxies occupy similar 
regions of each respective parameter space to the passive galaxies 
with only a few star forming outliers, typically with [Fe/H] ~ -1. 
The fitted slopes for [ C/Fe], [N/Fe] and [C a/Fe] are all in agreement 
with those found bv lGraves et al. I J2007h within the joint parame- 
ter uncertainties. Th ey are somewhat larger than those reported by 
ISmith etHll d2009ah for their simple two-parameter fits to Coma 
dwarf galaxies. However, those fits may well be poorly constrained 
given the narrow a range covered by their dwarf sample. 

4.5 Trends with Radius 

In Fig. [T3] we plot the measured stellar population parameters 
against cluster-centric radius for three out of four subsamples, only 
excluding the star forming systems again for clarity. The trend lines 



The Stellar Populations of Bright Coma Cluster Galaxies 15 




Figure 12. The six derived stellar population parameters as a function of galaxy velocity dispersion for the 222 passive and 94 emission line galaxies fit by 
Ez-AGES. Symbols are the same as in Fig. [5] The legend again contains the slope of a fit to the passive sample, first not masking and then masking the blue 
subset. The latter fit is also plotted as a red line in each panel. The blue line is a fit to the red AGN/LINER sample identified in the Section |431 Error bars are 
median parameter errors for the passive, red, and emission, blue, samples. The black dashed lines in the first two panels denoted the age (17.7 Gyr) and [Fe/H] 
(+0.2) limit of the Schiavon models. 




Figure 13. Stellar population parameters as a function of projected cluster-centric radius for the passive and emission line galaxies that were fit by Ez-AGES. 
Symbols are as in Fig. [6] Trend lines are fit to the raw data which are then binned with ~ 27, ~ 1 1 and ~ 7 galaxies per bin for red, cyan and black 
bins respectively. Star forming galaxies are not plotted for clarity. The legend contains the slope of the passive galaxy fit, excluding the blue subset, and its 
significance. 
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Figure 14. Stellar population parameter residuals as a function of cluster-centric radius. Symbols as in Fig.[6]and binning as in Fig. 1 1 3| Trend lines are fit to 
the raw data. The legend contains the slope of the passive galaxy fit, excluding the blue subset, and its significance. Star forming galaxies are not plotted for 
clarity. 



Table 5. Stellar population parameter gradients per degree, from the linear fits presented in Fig. ll3l and the planar fits displayed in Fig. 1141 The number of 
galaxies in each respective fit is displayed in parentheses in the table header and the significance of each fit in parentheses in the body of the table. For the 
passive sample the results in the first column include the blue subset while second column do not. 



Parameter 


Passive (222/220) 


AGN/LINER (44) 


SF (24) 


AED (26) 


Parameter - R cc (Fig. 1 1 3| 



log(age) 
[Fe/H] 

[Mg/Fe] 
[C/Fe] 
[N/Fe] 

[Ca/Fe] 



-0.074±0 
0.029±0 

-0.042±0 
0.023±0 

-0.061±0 

-0.014±0 



019(3.8) 
014(2.1) 
012(3.5) 
012(1.9) 
014(4.4) 
015(1.0) 



-0.082±0. 

0.029±0. 
-0.043±0. 

0.023±0. 
-0.064±0. 
-0.014±0. 



018(4.4) 
014(2.0) 
011(3.8) 
012(1.9) 
015(4.4) 
015(0.9) 



0.154±0.033(4.6) 
0.004±0.025(0.2) 
0.006±0.019(0.3) 
0.025±0.023(1.1) 
0.025±0.024(1.0) 

0.007±0.026(0.3) 

^Parameter - Rrr (Fig. 1 141 



0.079±0.067(1.2) 
0.097±0.061(1.6) 
0.061±0.056(1.1) 
0.069±0.055(1.3) 
-0.269±0.060(4.5) 
-0.116±0.076(1.5) 



-0.134±0.046(2.9) 
0.019±0.034(0.6) 
-0.026±0.028(0.9) 
-0.009±0.030(0.3) 
-0.029±0.034(0.9) 
-0.027±0.035(0.8) 



log(age) 
[Fe/H] 

[Mg/Fe] 
[C/Fe] 
[N/Fe] 

[Ca/Fe] 



-0.041±0 
0.029±0 

-0.033±0 
0.037±0 

-0.042±0 

-0.012±0 



018(2.3) 
015(1.9) 
011(2.9) 
012(3.1) 
014(3.0) 
015(0.8) 



-0.047±0. 

0.028±0 
-0.033±0. 

0.039±0. 
-0.044±0. 
-0.012±0. 



018(2.6) 
016(1.8) 
010(3.3) 
012(3.3) 
014(3.1) 
016(0.7) 



-0.081 ±0.036(2.2) 
-0.011±0.025(0.5) 
0.022±0.022(1.0) 
0.014±0.024(0.6) 
-0.004±0.026(0.2) 
0.004±0.027(0.2) 



0.085±0.079(1.1) 
0.065±0.061(1.1) 
0.076±0.057(1.3) 
0.070±0.052(1.3) 
-0.229±0.068(3.4) 
-0.130±0.079(1.6) 



-0.090±0.053(1.7) 
0.004±0.042(0.1) 

-0.004±0.034(0.1) 
0.001 ±0.033(0.0) 

-0.010±0.041(0.3) 

-0.019±0.049(0.4) 



Table 6. Stellar population parameter radial gradients obtained by fitting a plane to each passive galaxy subset as defined by the first column. Column 2 gives 
the number of galaxies in each fit. Significant gradients (> 2a) are denoted by bold text. ETG and LTG refer to early-type and late-type galaxies respectively. 





Criteria 


Ngal 


log(age) 


[Fe/H] 


[Mg/Fe] 


[C/Fe] 


[N/Fe] 


[Ca/Fe] 




All 


220 


-0.047±0.018 


0.028±0.016 


-0.033±0.010 


0.039±0.012 


-0.044±0.014 


-0.012±0.016 


a < 


125 km s" 1 


119 


-0.006±0.025 


-0.006±0.023 


-0.033±0.019 


0.032±0.019 


-0.032±0.021 


0.007±0.025 


a > 


125 km s" 1 


101 


-0.084±0.025 


0.067±0.019 


-0.036±0.014 


0.047±0.016 


-0.072±0.018 


-0.025±0.019 




ETG 


158 


-0.039±0.024 


0.039±0.020 


-0.032±0.016 


0.048±0.017 


-0.092±0.020 


-0.018±0.019 


ETG a 


< 125 km s" 1 


79 


0.046±0.038 


0.016±0.036 


-0.050±0.032 


0.069±0.029 


-0.105±0.035 


-0.005±0.042 


ETG a 


> 125 km s' 1 


79 


-0.092±0.033 


0.061±0.023 


-O.O35±0.018 


0.034±0.020 


-0.099±0.022 


-0.029±0.024 




LTG 


20 


-0.086±0.067 


-0.029±0.053 


-0.024±0.045 


-0.005±0.044 


0.011±0.050 


-0.022±0.055 
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are linear least-square fits to the raw data, weighted by their uncer- 
tainties, which are then binned for clarity. Each bin contains ~ 28, 
~ 1 1 and ~ 7 galaxies per bin for the passive, AGN/LINER and 
AED samples respectively. The results of these fits, including those 
to the star forming galaxies, are presented in Table[5] first section. 

For the quiescent galaxies significant correlations are obtained 
for log(age) (3.8a), [Mg/Fe] (3.5a) and [N/Fe] (4.4a). These trends 
typically get marginally stronger when masking the blue subset. No 
significant trend is seen for [Fe/H]. This implies passive galaxies 
on the outskirts of the cluster are, on average, 0.15 dex (~ 2 Gyr) 
younger and have more extended star formation histories than those 
in the cluster core. Both AGN/LINER galaxies and those with un- 
classified emission display significant age trends with comparable 
slopes. It is only the small sample of star forming galaxies which do 
not possess a correlation between SSP-equivalent age and radius. 

As demonstrated for the index-radius relations in Fig. |7J con- 
trolling for variations in velocity dispersion with radius can have a 
substantial effect on the significance of the observed radial trends. 
Therefore in Fig. [14] we plot the residual stellar population parame- 
ters, again derived by fitting a multivariate plane to the data, against 
radius and as such track the parameter variations at fixed velocity 
dispersion. The results of these fits are tabulated in the second sec- 
tion of Table[5] 

For the passive galaxies, significant, although typically 
weaker, trends are still found for log(age) (2.6<r), [Mg/Fe] (3.3a), 
[C/Fe] (3.3a) and [N/Fe] (3.1a). Thus bright quiescent galaxies on 
the outskirts of the cluster are on average ~ 20% younger with 
lower [Mg/Fe] and [N/Fe] but higher [C/Fe] than those in the core 
at fixed velocity dispersion. 

By contrast the only remaining significant trends for the emis- 
sion line galaxies are log(age) (2.2a) for AGN/LINER systems 
and [N/Fe] (3.4a) for star forming galaxies. This may well stem 
from the small size of these subsamples. In fact we find that if 
the AGN/LINER and AED galaxies are fit together, an approach 
that is not wholly unjustifiable given the identified similarities be- 
tween the two samples, a log(age) correlation of — 0.101±0.028 
dex deg -1 is recovered. Repeating this for [Fe/H] and [Mg/Fe] we 
find no significant trends for the combined sample. 

Again, a valid question here is to what extent are the rela- 
tively weak but significant radial trends obtained for the red pas- 
sive sample driven by the outermost data points. Once more we 
have repeated the fits with red passive galaxies beyond R cc = 1.55° 
masked. For the planar fits the results for [Mg/Fe] (2.5a), [C/Fe] 
(3.7a) and [N/Fe] (3.3<r) remain comparatively unchanged in terms 
of strength and sign. The same is also true for the linear fits, al- 
though [C/Fe] now becomes significant at the 3.4a level. The ra- 
dial trend in log(age), however, becomes insignificant in the lin- 
ear (1.7(f) and planar (1.2a-) fits yet the signs remain the same as 
those presented in Table [5] Thus it would seem the trend toward 
younger ages with radius begins to become significantly defined at 
Rcc > 1-5° for the red passive galaxies, although the present data 
cannot trace the form of the relation in any detail. 

Next we seek to check the correlations found for the passive 
galaxies against other possible methods of subdivision such as mor- 
phology and a. The former is of particular interest since, as no mor- 
phological criteria were imposed, the morphology-density relation 
may well play a key role in driving the observed radial stellar popu- 
lation gradients for this sample. Thus in Table|6]we apply a variety 
of further cuts to the quiescent galaxies excluding the two particu- 
larly blue early- type galaxies. 

Firstly, the sample is split into subsets based on their veloc- 
ity dispersion into low (a < 125 km s -1 ) and high (a > 125 km 



Table 7. Comparison between the results obtained for this work's passive 
galaxies, both full sample and core-South- West, and those in the literature. 
S06 gradients are again converted to per degree assuming R200 = 1.59° for 
Coma and we assume their [a/Fe] gradient is comparable to our [Mg/Fe] 
gradient. 



Source 


log(age) 


[Mg/Fe] 


This Work: All Passive 


-0.047±0.018 


-0.033±0.010 


This Work: Passive Core-SW 


-0.041±0.048 


-0.047±0.028 


S08 (Coma dwarfs) 


-0.369±0.069 


-0.128±0.039 


S06 


-0.045±0.011 


-0.030±0.006 



s _1 ) velocity dispersion bins. Here we see that, most notably, sig- 
nificant trends in log(age) (3.2a) and [Fe/H] (3.3a) are recovered 
for the high a sample. Next we separate the sample into early-type 
and late-type galaxies based on their NED morphologies. Here we 
find that the former does not have a significant radial age gradient, 
though it is consistent with that found for the passive galaxies as a 
whole taking into account the fact that these two samples are not 
independent. For the latter we find a sizable but insignificant nega- 
tive trend. At face value one might now be tempted to conclude that 
the gradient found for the entire sample (-0.047±0.018 dex deg -1 ) 
may well be driven by a combination of a very weak trend in the 
early-type galaxies modulated by a larger negative gradient in the 
late-type galaxies. To test this we subdivide the early-type galaxies 
into low and high a bins. Interestingly we find that the higher a 
early-type galaxies possess a strong negative age gradient, compa- 
rable in magnitude to the late-types, while the low a early-types 
show no significant gradient. This perhaps implies that the radial 
variation of the stellar populations of early-type galaxies in Coma 
may be more complex than the simple planar model applied here. 
Nevertheless we confirm that the morphology-density relation does 
not drive the log(age)-radius correlation seen for the passive sam- 
ple. 

In Table |7] we present a comparison of the significant results 
obtained in this work for the passive galaxies to those studies fea- 
tured previously in Table [3] excluding C02 who did not derive 
log(age) or [Mg/Fe] gradients. Again we fit separately to the Core- 
SW region subsample. We consider the [a/Fe] result of S06 com- 
parable to our [Mg/Fe] gradient. From this table it is immediately 
apparent that highly consistent age and abundance ratio gradients 
are found by this work and S06, despite each using different SSP 
models. It also indicates that S08 derived much stronger trends for 
their dwarf sample. We confirm that for bright Coma galaxies a 
limited azimuthal coverage does not result in significantly stronger 
gradients, again taking into account that the overall and Core-SW 
samples are not independent. 

4.6 The Z-plane 

It was lTrager et al . ( 2000b) who first highlighted an anti-correlation 
between age and metallicity at fixed velocity dispersion in early- 
type galaxies, terming it the Z-plane. In physical terms this rela- 
tion implies that younger galaxies, in an SSP-equivalent sense, have 
played host to a larger number of stellar generations allowing for 
greater metallic enrichment relative to older galaxies. To test for 
such a trend in our data we fit our quiescent sample with an Fe- 
plane of the form, 

[Fe/H] = a log a + b log age + c 
where the coefficients a and b give the dependence of [Fe/H] on a 
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Figure 15. Two edge on views of our version of the lTrager et alj {2000b) Z-plane, plotted here as an Fe-plane, for our quiescent sample. In the right panel we 
plot a typical 68% confidence error ellipse which, as discussed in the text, is nearly parallel to the observed age-metallicity anti-correlation. The dashed line in 
this panel shows the best fit relation. 
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Distance from mean value along best fit line 

Figure 16. Distribution of distances from the mean displacement following 
each point in the right hand panel of Fig.[l5]being deprojected onto the best 
fitting line. The one-dimensional \ 2 , degrees of freedom and q, the proba- 
bility that the observed distribution could be attributed solely to correlated 
errors, are shown. 

at fixed age and on age at fixed a respectively. The fit is weighted 
by errors in [Fe/H] and age but no attempt to account for corre- 
lated errors is made at this stage (but see below). We recover a = 
0.93±0.10, b = -0.90±0.08 and c = -1.24±0.17 with an rms scatter 
of 0.17 and in Fig.[l5]plot two projections of the fitte d plane. The 
slope o f [Fe/H] with a is stee per than that found by iTrager et al.l 
fcOOObh . a = 0.48±0. 12, and bv lSmith et at] J2009ah for their Coma 
dwarf sample, a = 0.51 ±0.2. However, the age-metallicity correla- 
tion we derive is comparable within the joint errors to that obtained 
by both works, b = -0.74±0.09 from the former and b = -0.74±0. 14 
from the latter. We estimate an intrinsic scatter of 0. 13 in [Fe/H] and 
in the left panel of Fig. [15] just by e ye, note an increased scatter for 
lower a galaxies, as also found bv lSmith et alj |2009a). Quantita- 
tively, the intrinsic scatter increases from 0.04 for log a > 2. 1 to 
0.17 for logo- < 2.1. 

Correlated errors between age and [Fe/H] can be a serious 



problem when trying to establish whether an age-met allicity anti- 
correla tion exists at fixed velocity dispersion (e.g. ITrager et al.1 
2000b). Indeed it may be the case, given the fact that the orientation 
of the error ellipse in the right hand panel of Fig. [15] is parallel to 
the slope of the dashed line, that the observed anti-correlation can 
be completely accounted for by uncertainties in the data. To test 
wheth er this is the case we employ a similar method to Rawle et al. 
( 2008) and deproject the data onto the fitted plane (the trend line 
in Fig. Q3] second panel). Next we compute the one-dimensional 
X 2 for the data along the line about a constant value, taking the 
one-dimensional error to be the semi-major axis of the la (40% 
joint confidence in two-dimensions, 68% in one-dimension along 
the major axis) error ellipse obtained from our error model simu- 
lations rescaled to the S/N of each data point. See Appendix |B1 for 
an examination of the robustness of this choice. In Fig. [16] we plot 
the distribution of one-dimensional residuals against an arbitrary 
y-axis, together with the one-dimensional x 2 > degrees of freedom 
and Q, the probability that the observed distribution could be at- 
tributed solely to correlated errors. This test demonstrates that the 
Z-plane has a high probability of being resolved and would require 
an increase in errors in log(age) and [Fe/H] of >20% to become 
unresolved, which, given the analysis conducted in Appendix [B] is 
very unlikely. 



It is also possible that the correlated errors may bias 6, the 
slope of the age-metallicity trend at fixed a. To test for this effect 
we use simple simulations that take as their input the observed ages 
and velocity dispersions and new artificial [Fe/H] values which are 
assigned using the Z-plane with coefficients as stated above. We 
generate new realisations of this dataset by resampling the ages 
and metallicities using the error ellipse, again rescaled to the S/N 
of each data point. To each new realisation we fit the Z-plane and 
obtain a new b coefficient. Repeating this process 500 times we find 
that, for an anti-correlation as strong as the one we observe, b is not 
significantly influenced by the correlated errors between age and 
metallicity. 
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5 DISCUSSION 

5.1 Downsizing in the Coma Cluster 

Evidence of downsizing, referred to here as a correlation between 
galaxy velocity dispersion and SSP-equivalent age, has been iden- 
tified by many recent studies on t he stellar populations of quiescent 
galax i es in the loca l unive r se (e.g. Caldwell et al. 2003; Nelan et al. 
20051: ISmith et all 120061: iGraves et all |2007| ; ISmith et alj 12007k 
Thomas et a l. 2009; G raves et alj|2009l) . This relation can be con- 



sidered the archaeologic al equivalent of the downsizing trend ob- 
served by authors such as |Cowie etaljjl996h and more recently by 
IPerez- Gonzalez et al. ( 2008), who find that the total luminosity of 
galaxies which are undergoing rapid star formation declines as a 
function of redshift. 

The slope of the age-cr relation reported by those works 
listed above are between 0.35 dGraves et al]|2007l) and 0.72 (S06) 
across approximately the same cr range as studied here. In fact 
iBernardi et alj |2006) find a slope of 1 . 1 based on stacked early- 
type galaxy spectra from the SDSS. This sizable spread can, at 
some level, likely be attributed to a host of systematics such as sam- 
ple selection (morphology, colour or both), data processing tech- 
niques and the SSP model of choice for each study. Nevertheless a 
picture is developing such that age does correlate with velocity dis- 
persion in the sense that higher cr quiescent galaxies formed their 
stars earlier tha n lower cr galaxies. 

Previously Sanchez-Blazqu ez et alj J2006h and T08 have re- 
ported that early-type galaxies in the core of the cluster show 
no signal of do wnsizing. T08 go o n to highlight that the Coma 
sample used inlThomas et aU fc005h . 35 early-type galaxies from 
Mehl ertet alj d2000h . also does not support downsizing and that 
the ag e-cr relation t hey fin d for the 97 Coma red-sequence galax- 
ies in iNelan et alj J2005l . hereafter NFPS) is the result of no 
Bal mer line emission infill correction being applied. On the con- 
trary iMatkovic et al. I J2009h do find, albeit limited, evidence for 
downsizing in their Coma galaxies. Given the wealth of studies on 
the quiescent galaxy population in the local volume that support a 
significant age-cr relation, one must ask why evidence for its signa- 
ture in Coma seems inconclusive at best. 

With this in mind, in this section we seek to test the robust- 
ness of our observed age-cr correlation against a number of sys- 
tematic effects. These are the level at which we claim a galaxy to 
be quiescent, the approach used to correct the observed indices for 
each galaxy's velocity dispersion and the choice of SSP model em- 
ployed. 

As detailed earlier, our analysis requires that a galaxy must not 
have emission in any of H/3, [OIIIJ5007, Ha or [NII]6585 detected 
by GANDALF with a A/N ^ 4 in order to be considered quiescent. 
This detection threshold was set in accordance with the authors of 
the routine but we note that T08 opt for a much stricter cut, clean- 
ing any emission observed at the A/N > 2 level. In part this more 
severe criteria likely helps compensate for their study not having 
access to the, typically, stronger Ha and [NII]6585 lines. That said 
the S/N of their spectra, which is much higher than that used here, 
will result in more sensitivity to emission and so allow the detection 
of fainter lines at fixed A/N. To test what effect cutting at this level, 
which amounts to a detection threshold of EW(H/3) ~ 0.13 A in 
our SDSS spectra, has on our observed downsizing trend we re-fit 
the 158 of the 222 passive galaxies that meet this criterion and have 
measured stellar population parameters. We find age oc a °- 96±0 - ' > 
consistent with the slope obtained when fitting the entire sample. 
Thus a stricter emission line cut appears not to affect the downsiz- 
ing trend we recover. 



While perhaps only a minor effect, iKelson et alj J2006h find 
some evidence for the particular scheme used to velocity broaden- 
ing correct their index data influencing their measured parameter 
correlations. Here we opt to test two further correction techniques, 
that emplo yed by LlCK_Ew a nd the more standard multiplicative 
method of Trag er et al. Ul998h . The former smooths the observed 

spectra by a smooth = J af ick - a? es - a 2 gal where au ck is the 
Lick resolu tion of the partic ular index being measured (taken from 
Table 1 of ISchiavonl 12007). a res is the SDSS instrumental reso- 
lution and <jgai the galaxy velocity dispersion. Should the sign of 
the term under the square root turn negative the code u ses multi- 
plicati ve corrections based on smoothed SSP spectra from lSchiavonl 
(2007) . We take the veloc ity broadening polynomial coefficients 
for the lTrager et alj j 19981) multiplicative method from their Table 
5. The technique used by LlCK_Ew yields age oc a °- S3±0 ' a7 while 
the multiplicative route gives age oc (j - 7>i±0 - 0(i _ Thus, while the 
latter is a shallower trend, our correction scheme and both of those 
tested produce consistent results within the joint errors. 

To test the model dependence of our observed downsizing 
trend we make use of the empirically derived transformations of 
Smit h et al. (2009a) to convert our stellar population parameters to 
the TMBK model system. Fitting the age-cr relation to the trans- 

n -\- n oft 

formed data we find age oc a ' . Again this is consistent, 

within the joint uncertainties, to the result obtained when the Schi- 
avon models are employed. 

In summary, following the above tests, and those conducted 
earlier, we conclude that our SDSS data does robustly support 
downsizing in both early-type and passive galaxy populations in 
the Coma cluster. As demonstrated in Appendix lAl we have found 
good consistency between this work and NFPS in terms of derived 
age-cr relations. We speculate that small sample size may well ex- 
plain the absence of downsizing seen by T08 using their Keck/LRIS 
data, given the sizable scatter in the first panel of Fig.[T2] Further- 
mo re we have identified that, for galaxies which are in both studies, 
the lMoore et alj d2002l . hereafter M02) dataset results in a signifi- 
cantly flatter H/3-cr slope than this work. This issue, with as yet 
unknown origin, may well drive T08 to conclude the M02 data dis- 
plays no age-cr correlation. 



5.2 Environment 

The Coma cluster is an ideal place to study how environment af- 
fects the stellar populations of galaxies across a wide range of local 
densities. Its core represents one of the highest density regions that 
is readily observable for line strength studies while the cluster's 
outskirts provide excellent insight into galaxies that are, perhaps 
for the first time, just starting to feel the effects of a cluster envi- 
ronment. Indeed evidence has been mounting for some time that 
galaxies in the cluster's core display systematically different line 
indices, and by inference stellar populations, than those in the outer 
regions. 

The first study to highlight the spatial dependence of the stel- 
lar population s of C oma members using line indices was that of 
iGuzman et ail |l992). They observed 51 Coma elliptical galaxies 
located in the cluster's core and "halo", defined by them having 
a projected radius of 1° < R < 5°, with cr > 100 km s -1 . As 
such their sample is comparable to the high cr end of this work. 
They found that their Mg2 -cr relation displayed a small environ- 
mental dependence such that, at fixed cr, galaxies in the core have 
higher index strengths than those in the cluster's halo. Employing 
early stellar population models, they went on to attribute this trend 
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to variations in age with younger galaxies at larger cluster-centric 
radii. 

Sh ortly after the work o f Guzman et al., the aforementioned 
study o flCaldwell et al ] dl993h uncovered a sample of post-starburst 
early-type galaxies in the outer South- West region of the cluster. 
These galaxies were identified based on their strong high-order 
B aimer lines, indicative of a subpopulation of A stars with lifetimes 
~ 1 Gyr. Caldwell et al. could not comment further on whether 
this post-starburst population were localised to the NGC 4839 sub- 
group or spatially a more general phenomenon, as they were only 
able to observe one fi eld on the clusters outskirts (although see 
ICaldwell & Rose! 1 19971) . That said, very recently M ahaian et al.l 
d2010h found that k+a galaxies in Coma are in fact uniformly dis- 
tributed about the cluster, only avoiding regions with substantial 
X-ray emission such as the core. 

C02 once more addressed the effect of environment on the 
stellar populations of Coma members, this time spanning a wide 
range in total magnitude from dwarf to giants, using Lick indices. 
They found evidence for a significant radial gradient in Mg2, even 
after factoring out the index's correlation with luminosity. As pre- 
viously mentioned, based on the weaker trends they recovered for 
HP, they interpreted this as a metallicity gradient. More recently, 
S08 have concluded that the radial index trends in a sample of red- 
sequence Coma dwarf galaxies in fact support a strong gradient in 
their SSP-equivalent ages, with older dwarfs preferentially found in 
the cluster's core. Again neither of these studies were able to com- 
ment on the generality of this gradient due to limited azimuthal 
coverage. 

The results for our quiescent sample, presented in Table [5] 
provide further support for an interpretation of the observed radial 
index trends in terms of variations in galaxy age rather than metal- 
licity. We also find evidence for significant gradients in [Mg/Fe], 
[C/Fe] and [N/Fe]. These findings indicate that passive galaxies on 
the outskirts of the cluster are younger, formed their stars over a 
longer period of time and have higher [C/Fe] but lower [N/Fe] than 
galaxies located in the cluster's core. Furthermore, when consider- 
ing our core-South- West subsample we find no significant increase 
in recovered age or [Mg/Fe] gradients, providing the first hint that, 
at least for bright galaxies, the South- West field is not unique. 

In Tables[3]and[7]we have shown that the samples of S08 and 
C02 typically display much stronger cluster-centric radial gradients 
in terms of line indices and stellar population parameters than de- 
rived for our sample. We have proposed that this finding can be 
attributed to the differing luminosity coverage of the three samples 
rather than a bias being introduced by S08 and C02 only observing 
the cluster's South- West region. Taken together this implies that ei- 
ther; brighter galaxies had finished forming the majority of their 
stars before entering the cluster, star formation within these sys- 
tems was less affected by the cluster environment, or some mix 
of the two. The former relies on such ga laxies having undergone 
substantial pre-pr ocessing in groups (e.g iKodama & Smailll200fl ; 
iMcGee et alj|2009t) prior to being accreted into the cluster while 
the latter can likely be put down to their deeper potential wells 
being better able to retain gas for star formation. Irrespective of 
the precise scenario, our results indicate that the stellar populations 
of bright galaxies are less influenced by entry into the cluster en- 
vironment than those of dwarf galaxies. Furth e rmore , this picture 
fits well with that discussed by iGavazzi et al.l q201Cr> . They stud- 
ied galaxies in the Coma supercluster and report that the colours 
and luminosity distributions of high luminosity galaxies display 
little environmental dependence relative to those of low luminos- 
ity galaxies. Thus they propose a scenario where the former were 



in place at earlier times while the latter are currently, cosmologi- 
cally speaking, undergoing environmentally driven transformations 
as they enter denser structures. 

We have also compared our findings to those from S06 who 
studied quiescent galaxies in a host of nearby clusters across a 
comparable a range to this work. Using a similar planar analysis 
technique to that used here but employing different SSP models, 
they found average cluster-centric radial gradients that are highly 

consistent with those reported here. 

In ter ms of other recent mo re general studies. IBernardi et al.l 

(2006) and lClemens et al.l {2009) both used large samples of early- 
type galaxies selected from the SDSS to test the environmental de- 
pendence of their stellar populations. Using differing techniques, 
they split galaxies into high and low local galaxy density subsam- 
ples. Here we note that the outskirts of Coma are unlikely to be 
comparable to low density or field regions due to the pre-processing 
of in-falling galaxies. Nevertheless, both studies report that early- 
types in low density regions are systematically younger than those 
in high density environments, in agreement with our results. By 
contrast, Idemens et al. I d2009h find marginal evidence that [Z/H] 
is lower in den se regions but that [q /Fe] does not depend on envi- 
ronment, while IBernardi et al.l ( 120061) identify no [Z/H] correlation 
but conclude that a-enhancement is greater in high density regions. 
Thus it would appear the latter study fits better with our findings. 

Finally, we are able to compare our results with the predictions 
from t he semi-analytic galaxy formation models of lDe Lucia et al.l 
(2006). Their Fig. 8a indicates how galaxy median luminosity- 
weighted age varies, on average, as a function of physical distance 
out from cluster centres. While direct comparison is complicated 
by deprojection effects, mass segregation and conversion from 
lumin osity- weighted to SSP-equivalent age 1 Trager & Somervilld 
2009), this figure predicts a typical change of ~ 13% between 
O.IR200 and R2oo- This trend is consistent with our finding of 19% 
± 8% within the same radial range in Coma. In addition, our re- 
sults quali t ativel y agree with those from the N-body simulations of 
iGao et all {2004). They find that subhaloes located near the centre 
of cluster-sized parent haloes were accreted earlier. Thus a combi- 
nation of this correlation with any mechanism that acts to inhibit 
star formation in galaxies as they enter the cluster environment will 
lead to a cluster-centric radial gradient in their stellar populations. 



6 CONCLUSIONS 

In this paper we have studied the stellar populations of bright Coma 
cluster galaxies using SDSS archival spectroscopy and up-to-date 
SSP models. The key findings of this work can be summarised as 
follows: 

(i) Of the 110 galaxies (31% of the entire sample) with de- 
tectable emission we find 46% show AGN/LINER like emission, 
29% have ongoing star formation and 25% have an ambiguous 
emission driver based on this data. The fraction of AGN/LINER 
systems is seen to increase with increasing galaxy luminosity while 
the star forming fraction decreases. 

(ii) For the 246 passive galaxies that meet our selection criteria 
we find strong correlations between absorption line index strength 
and velocity dispersion for CN2, C4668, Mgb and H/3. Particularly 
in the case of H/3, scatter about these relations is seen to increase 
with decreasing velocity dispersion. At low a the emission line 
galaxies typically have stronger H/3 and weaker metal lines than 
quiescent galaxies, while at high a both subsets have comparable 
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index strengths. This trend is driven by the aforementioned tran- 
sition from star forming to AGN/LINER galaxies with increasing 
luminosity, and therefore a. 

(iii) Using a planar analysis technique that factors out index cor- 
relations with a, we find significant cluster-centric radial gradients 
in HyS, Mgb and C4668 for passive galaxies. Comparing our re- 
sults with those of other recent studies that typically show stronger 
radial index gradients, we have shown that this difference can be 
attributed to the luminosity range covered by each study and not 
the uniqueness of the NGC 4839 in-fall region. 

(iv) We have demonstrated that the vast majority of our passive 
sample do not play host to a frosting of young (< 1 Gyr) stars on 
top of an underlying, mass dominant old population. 

(v) For 222 passive galaxies in Coma we find a strong positive 
correlation between their SSP-equivalent age and velocity disper- 
sion. The scatter about this relation increases toward lower a. This 
trend is found to be robust against variations in sample selection 
criteria (early-type vs quiescent), emission line detection thresh- 
olds, index velocity broadening corrections and the particular SSP 
model employed. Weaker positive correlations are recovered be- 
tween a and [Fe/H], [Mg/Fe], [C/Fe] and [N/Fe]. Thus higher a 
quiescent galaxies in Coma are found to have formed their stars ear- 
lier, more rapidly and with higher [Fe/H], [C/Fe] and [N/Fe] than 
cluster members with lower velocity dispersions. 

(vi) Galaxies which show AGN/LINER emission are typically 
found to have younger SSP-equivalent ages than quiescent galaxies 
at fixed velocity dispersion. 

(vii) Significant cluster-centric radial stellar population gradi- 
ents are found for our passive sample in terms of their SSP- 
equivalent age, [Mg/Fe], [C/Fe] and [N/Fe]. These trends are in the 
sense that, at fixed a, galaxies on the cluster's outskirts are younger, 
formed stars over a longer period of time and have higher [C/Fe] 
but lower [N/Fe], Again we attribute fainter luminosity coverage 
rather than the uniqueness of the NGC 4839 group as the driver for 
the stronger radial gradients found by other recent studies. We find 
that the morphology-density relation does not drive the recovered 
age-radius relation. 

(viii) Interpreting galaxies with barely detectable emission as 
weak AGN/LINER systems, galaxies with this emission driver 
are found to have a significant cluster-centric radial gradient in 
their SSP-equivalent ages. This trend is more than a factor of two 
stronger than that seen for quiescent galaxies. 

(ix) Finally, for our passive sample we have found an age- 
metallicity anti-correlation that is consistent with previous works. 
We have shown that this anti-correlation is resolved, and so is not 
the result of correlated errors between age and metalicity, and that 
correlated errors do not significantly bias the measured slope. 

We have demonstrated that the stellar populations of bright 
quiescent galaxies in Coma depend strongly on their a and, rela- 
tive to dwarf cluster members in the South- West, less so on their 
local environment. We have also shown that the stellar populations 
of such galaxies in the NGC 4839 in-fall region are not unique 
with respect to galaxies elsewhere in the outskirts of the cluster. 
In a forthcoming paper, we will explore these two related issues 
using an enlarged sample of faint galaxies. Firstly, to determine 
whether the strong environmental trends seen for dwarf galaxies 
in the South- West are localised there or a more general effect in 
the outskirts of the cluster. Secondly, to identify the transition a, 
or range of velocity dispersions, at which the stellar populations 
of a galaxy start to become significantly more affected by its local 
environment. 
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APPENDIX A: INDEX-INDEX COMPARISONS 

In this section we compare our measured index data with other 
studies from the literature both focusing on and involving Coma 
cluster galaxies. Specifically we opt for two datasets which possess 
a sizable number of galaxies in common with this work, namely 
M02 and NFPS. The former has 78 galaxies in common with our 
passive sample while the latter has 75, permitting a suitable ro- 
bustness check of our measurements. Here we make no attempt to 
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match spectroscopic fibre sizes between datasets, noting that M02 
observed with a 2.7" diameter fibre and NFPS with a 2" diameter 
fibre. In addition we note that both studies are based on flux cali- 
brated spectra without explicit correction to the Lick/IDS system, 
in a similar vein to our data. 

In Fig. I A 1 1 and Fig. IA2I we present diagrams comparing the 
six indices used in our analysis with those from the literature stud- 
ies. For the M02 sample not all indices are available due to the 
wavelength coverage of that work's spectra. Mean index offsets are 
computed as <Literature - This work>. The quoted reduced \ 2 is 
computed for a model with no offset if the offsets significance is 
< 3a, otherwise it is about a model taking into account the mean 
offset. A one to one compatibility between the index data is also 
represented by the straight line in each panel. 

From Fig. lAll we see that only C4668 has a significant offset at 
the 4<r level. H/3 and Mgb have relatively high reduced x 2 , indicat- 
ing perhaps underestimated errors in this work or M02. In Fig. lA2l 
significant offsets are seen for Ca4227 (6a) and C4668 (5a). When 
these systematic shifts are accounted for, however, both Ca4227 
and C4668 from NFPS and this work are in agreement. High re- 
duced x 2 is om y obtained for CN2, again probably indicating un- 
derestimated errors. 

To summarise, generally speaking good agreement is seen be- 
tween this work and the two literature overlap samples. There is, 
however, reasonable evidence for a systematic offset in our C4668 
index measurements relative to those in both M02 and NFPS. 

Next we seek to identify what is driving the apparent lack of 
downsizing in the M02 data, as determined by T08, and a weaker 
signal of its presence in the NFPS data relative to this work. The 
logical place to search for the source of these differences is H/3, the 
main SSP-equivalent age indicator in all three studies, as a function 
of velocity dispersion. With this in mind, in Fig. |A3| and Fig. lA4l we 
plot Aindex-cr for both M02 and NFPS where Aindex is computed 
as Literature index - This work and the velocity dispersion mea- 
surements are those from this work. 

Fig. IA3I shows significant Aindex-cr correlations for H/3 and 
C4668. Importantly this demonstrates M02 measure stronger H/3 at 
high a and weaker H/3 at low a than our SDSS data. This in turn im- 
plies, everything else being equal, they would recover younger ages 
for higher a galaxies and older ages for lower a galaxies than this 
work. Thus the age-a relation would become flattened, although it 
is unclear to what extent. By contrast, in Fig. IA4I we see that the 
NFPS overlap displays no similarly significant AH/3-ct correlation. 
Here we do note that, while consistent with zero, the NFPS slope 
is also consistent with the trend seen for the M02 sample. In this 
case we are able to take the comparison a step further and compare 
both H/3-a and, using the Schiavon models, log(age)-a relations. 
For the former we obtain the slopes H/3 n / ps -cr = -0.93 ± 0.09 and 
HPsdss-a = -0.99 ± 0.16, highly consistent within the uncertain- 
ties. For the latter we are restricted to the 58 galaxies in common 
that Ez-AGES is able to fit and recover slopes of log(age)jvFPS- 
a = 0.67 ± 0.13 and log(age) sdss -(j = 0.77 ± 0.14. Again these 
trends are consistent within the uncertainties. 

In summary, we have demonstrated that NFPS and our SDSS 
data imply compatible downsizing trends where they overlap. How- 
ever, for the M02 and SDSS overlap we find a significantly flatter 
H/3-cr relation from the former, which in turn may result in T08 
obtaining no age-cr correlation when employing the M02 dataset. 
Given the fact that we estimate velocity broadening differences in- 
troduce ~ 15% difference in H/3-a slopes, it is unclear what drives 
this discrepancy. 



Table Bl. Relevant parameters for the galaxies used to test the error model. 
The first two columns give the galaxy ID, for comparison with Fig. lB2l and 
signal-to-noise. The final two columns show the correlation coefficient (p) 
between log(age) and [Fe/H] and the orientation of the resultant error el- 
lipse (6), with respect to the log(age) axis, as computed from each galaxy's 
simulations. Galaxies el-4 were used to construct the error model. 
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Figure B2. H/3-<Fe> model grid with the galaxies used for the error 
model test highlighted by large black stars. For further details on the grid 
see Fig. II II 



APPENDIX B: ERROR MODEL ROBUSTNESS 

In this section we seek to test the robustness of the error model used 
to compute the stellar population parameter uncertainties both in 
terms of its ability to reproduce parameter error at a given S/N and 
suitably measure the correlated error between age and metallicity 
required for Section |4~6l To that end we select six galaxies of vary- 
ing S/N and at different positions on the H/3-<Fe> grid (see Fig. 
lB2tfrom our passive sample and generate 200 realisations of their 
index data using their respective uncertainties. Each realisation is 
ran through Ez- AGES using the same approach as described for the 
real data. Next we use the S/N = 55 (e4, see Table IB it simulation 
data, used to make the error model, to assess the age-metallicity 
covariance and construct the resultant error ellipse, which is in turn 
rescaled to the respective S/N of each of the six selected galaxies. 
Finally, we draw 1000 samples from the rescaled error ellipse and 
in Fig. lBll compare the predicted and actual age-metallicity distri- 
butions for the six galaxies. We also include the S/N = 25, 35 and 
45 galaxies used to build the error model since it is only the co- 
variance of the S/N = 55 galaxy which is currently being tested. In 
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Figure Al. Index-index comparisons between this work and data from lMoore et all (2002). No aperture corrections are applied and the data is compared at the 
Lick resolution. The mean offset is computed as <Literature - This work>. The quoted reduced \ 2 is about the equality line displayed in each panel unless 
the offset has a significance ^ 3cr, in which case the model takes into account the mean offset. For clarity we restrict the axis range of the Hbeta and Mgb 
panels, excluding one k+a galaxy in common between the two samples. Note this galaxy is included in the offset and x 2 calculations. 




Hbeta <Fe> Mgb 




CN2 Ca4227 C4668 



Figure A2. Index-index comparisons between this work and data from Nelan et al. ( 2005). No aperture corrections are applied and the data is compared at the 
Lick resolution. The mean offset is computed as <Literature - This work>. The quoted reduced \ 2 is about the equality line displayed in each panel unless 
the offset has a significance ^ 3a, in which case the model takes into account the mean offset. 
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Figure A3. Aindex-cr trends for the M02 overlap sample. Aindex is computed as Literature index - This work and the velocity dispersion measurements are 
those from this work. 
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Figure A4. Aindex-cr trends for the NFPS overlap sample. Aindex is computed as Literature index - This work and the velocity dispersion measurements are 
those from this work. 
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Figure Bl. Diagrams comparing the result of drawing 1000 samples from the rescaled error ellipse of e4, see Table IbTI (grey points) to the output of 200 
realisations of each galaxy's index data which has been ran through Ez-AGES (black points). As such this figure may be used to visually assess the suitability 
of the error model. The legend of each panel displays the galaxy ID from Table lBll and its S/N. 



addition in Table IB ll we present the relevant parameters for the six 
randomly selected galaxies and four error model galaxies. 

Fig. IB II demonstrates that the rescaled error ellipse does a 
good job at reproducing the scatter from the realisations. Quanti- 
tatively the error model is confirmed to better than 10% for these 
galaxies. Table lBll shows that p is stable to within 10-15% (< p > 
= -0.62 ± 0.06) and 6 varies by at most 30% « > = -34.42 ± 
5.35°). Indeed, by comparing Table lBll and Fig. lB2l we see that for 
the H/3-<Fe> parameter space occupied by the bulk of our sample 
galaxies 9 is between -24° and -37° . Given its associated p and 6 
values, which are close to the means computed above, we choose 
to use the S/N = 45 (el) as the base error ellipse for rescaling when 
testing the impact of correlated errors in Section l4~6l 



APPENDIX C: RADIAL TRENDS: TESTING A NULL 
MODEL 

Given the significant, but relatively weak, radial trends detected for 
the passive sample in terms of their indices and stellar population 
parameters it is interesting to test the validity of these results, and 
indeed the fitting procedure, against a null model. As such, in this 
section we undertake simulations based on an artifical version of 



this sample, excluding the blue subset, in which there is no radial 
correlation. As input data for this test we use the velocity disper- 
sions and half-light radii as measured and generate fake indices and 
stellar population parameters using the coefficients derived from 
our planar fits but with the radial component set equal to zero. 
Next, for each of the six indices and stellar population parameters 
we create 1000 realisations of the input data using the measure- 
ment errors and intrinsic scatter derived from the planar fits to the 
real data. Finally, each realisation is ran through the same planar 
fitting procedure used for the real dataset and in Fig. IC II and Fig. 
|C2| we histogram the output of these simulations overlayed with the 
position of the coefficient obtained from the real data. 



Fig. |C II and Fig. |C2| clearly demonstrate the ability of the fit- 
ting procedure to account for intrinsic scatter and, on average, re- 
turn unbiased coefficients. The standard deviations of the distribu- 
tions presented in both figures are slightly larger than the errors 
given in Tables [2] and [5] because here the intrinsic scatter is resam- 
pled as well as the measurement uncertainties. Nevertheless, the 
signifiance of the C4668, H/3 and Mgb index trends and log(age), 
[Mg/Fe], [C/Fe] and [N/Fe] stellar population parameter trends is 
apparent. 
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Figure CI. Distribution of radial coefficients (o in equation [3 derived for the null model in which the indices of the passive sample are not correlated with 
projected cluster-centric radius. The black vertical bar denotes the value of the coefficient measured from the real data. The fraction of trials that yield an \a\ 
larger than the observed value is noted in parentheses. 
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Figure C2. Distribution of radial coefficients (a in equation pi] but with / set equal to each stellar population parameter in turn) derived for the null model in 
which the stellar population parameters of the passive sample are not correlated with projected cluster-centric radius. The black vertical bar denotes the value 
of the coefficient measured from the real data. The fraction of trials that yield an \a\ larger than the observed value is noted in parentheses. 
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Table Dl. Absorption-line index and velocity dispersion measurements used in this work. Galaxies are identified by their position on the sky (J2000) and, 
where available, their GMP number. The quoted S/N is per A and the velocity dispersion measurements are not aperture corrected. The indices and their errors 
are given in conventional units, i.e. magnitudes for CN1, CN2, Mgl and Mg2, angstroms for all others. The flag column is binary with indicating no detected 
emission with A/N = 4 and 1 denoted galaxies with detected emission. The full version of this table will be provided in the electronic version of the journal. 



RA 


DEC 


GMP 


S/N 


a 


(Terr 


H<5A 


H<5F 


CN1 


CN2 


Ca4227 


G4300 


H7A 


H7F 


Fe4383 


Ca4455 


Fe4531 


C4668 


H/3 


Fe5015 


Mgl 


Mg2 


Mgb 


Fe5270 


Fe5335 


Fe5406 


Fe5709 


Fe5782 


H (5" .A. 




1 <o y ip 




Ca4227 err 


G4300 err 




H7F e rr 


Fe4383 err 


Ca4455 err 


Fe4531 err 


C4668 err 


H/3 err 


Fe5015 err 


Mglerr 


Mg2err 


Mgb err 


Fe5270 err 


Fe5335 err 


Fe5406 err 


Fe5709 err 


Fe5782 err 


Flag 










13:00:08.13 


+27:58:37.0 


2921 


46.7 


386.9 


6.9 


-2.63 


-0.04 


0.136 


0.179 


1.14 


5.40 


-6.83 


-2.04 


5.01 


1.68 


3.84 


7.91 


1.44 


5.30 


0.169 


0.325 


5.19 


3.06 


2.96 


1.97 


0.96 


0.95 


0.21 


0.14 


0.006 


0.007 


0.11 


0.18 


0.22 


0.14 


0.28 


0.16 


0.24 


0.36 


0.16 


0.34 


0.004 


0.004 


0.16 


0.19 


0.23 


0.18 


0.13 


0.11 













13:04:58.38 


+29:07:20.1 


282 


43.1 


299.7 


6.4 


-1.34 


0.64 


0.087 


0.112 


1.03 


5.08 


-5.85 


-1.64 


5.00 


1.30 


3.68 


7.90 


1.71 


4.77 


0.144 


0.288 


4.68 


2.76 


2.74 


2.00 


1.06 


0.76 


0.25 


0.17 


0.007 


0.008 


0.13 


0.21 


0.24 


0.15 


0.31 


0.17 


0.26 


0.39 


0.17 


0.37 


0.004 


0.005 


0.18 


0.21 


0.25 


0.19 


0.15 


0.12 













12:57:24.36 


+27:29:52.1 


4928 


48.3 


278.4 


4.5 


-1.93 


0.47 


0.103 


0.141 


1.13 


5.53 


-6.27 


-1.78 


5.16 


1.38 


3.58 


7.86 


1.28 


5.48 


0.147 


0.299 


4.82 


2.88 


2.70 


1.93 


0.97 


0.76 


0.22 


0.15 


0.006 


0.007 


0.12 


0.19 


0.23 


0.14 


0.29 


0.16 


0.23 


0.35 


0.20 


0.34 


0.003 


0.004 


0.16 


0.19 


0.22 


0.16 


0.12 


0.10 













13:01:53.73 


+27:37:28.0 


1750 


48.6 


277.1 


4.4 


-2.20 


0.19 


0.124 


0.161 


1.19 


5.48 


-6.80 


-2.00 


4.93 


1.47 


3.43 


8.01 


1.50 


5.55 


0.158 


0.316 


4.89 


3.19 


3.07 


2.05 


0.91 


0.94 


0.19 


0.13 


0.005 


0.006 


0.10 


0.17 


0.21 


0.13 


0.26 


0.15 


0.22 


0.34 


0.15 


0.33 


0.003 


0.004 


0.15 


0.19 


0.22 


0.17 


0.12 


0.10 













13:01:57.57 


+28:00:21.0 


1715 


40.9 


277.0 


6.1 


-2.71 


-0.23 


0.134 


0.171 


1.03 


5.05 


-6.09 


-1.92 


4.47 


1.04 


3.44 


8.44 


1.83 


4.80 


0.165 


0.316 


4.96 


3.13 


2.95 


1.71 


0.77 


0.99 


0.27 


0.18 


0.007 


0.009 


0.14 


0.22 


0.26 


0.16 


0.34 


0.18 


0.27 


0.40 


0.18 


0.39 


0.004 


0.005 


0.19 


0.22 


0.26 


0.20 


0.17 


0.14 


1 










12:59:35.71 


+27:57:33.3 


3329 


43.7 


273.5 


5.1 


-2.72 


0.14 


0.125 


0.166 


1.23 


5.49 


-6.24 


-1.86 


5.00 


1.36 


3.65 


6.87 


1.54 


5.44 


0.147 


0.292 


4.80 


2.99 


2.81 


1.95 


0.84 


0.65 


0.26 


0.18 


0.007 


0.009 


0.13 


0.22 


0.26 


0.16 


0.32 


0.17 


0.26 


0.39 


0.17 


0.37 


0.004 


0.005 


0.17 


0.20 


0.24 


0.18 


0.15 


0.13 













12:59:03.90 


+28:07:25.3 


3792 


46.0 


273.2 


4.8 


-2.32 


0.14 


0.137 


0.175 


1.13 


5.58 


-6.83 


-2.19 


5.23 


1.35 


3.39 


8.33 


1.23 


5.83 


0.163 


0.327 


5.20 


3.06 


2.77 


1.86 


0.86 


0.84 


0.22 


0.15 


0.006 


0.007 


0.11 


0.19 


0.23 


0.14 


0.28 


0.16 


0.24 


0.36 


0.16 


0.35 


0.004 


0.004 


0.16 


0.19 


0.23 


0.18 


0.13 


0.11 













13:01:33.60 


+29:07:50.1 


1990 


40.3 


269.6 


5.9 


-2.25 


0.19 


0.113 


0.145 


1.06 


5.80 


-6.36 


-1.84 


4.56 


1.52 


3.57 


6.01 


1.66 


4.56 


0.132 


0.279 


4.76 


2.75 


2.72 


1.87 


0.82 


0.72 


0.29 


0.19 


0.008 


0.009 


0.15 


0.24 


0.28 


0.18 


0.36 


0.19 


0.28 


0.43 


0.18 


0.40 


0.004 


0.005 


0.19 


0.22 


0.26 


0.20 


0.15 


0.13 


1 










12:56:18.60 


+26:21:32.0 




48.2 


263.2 


4.2 


-2.15 


0.16 


0.111 


0.151 


1.15 


5.16 


-6.12 


-1.70 


4.83 


1.32 


3.58 


8.70 


1.37 


5.33 


0.151 


0.298 


4.98 


3.09 


2.64 


2.01 


0.89 


0.90 


0.24 


0.16 


0.006 


0.008 


0.12 


0.20 


0.23 


0.14 


0.29 


0.16 


0.24 


0.34 


0.15 


0.33 


0.004 


0.004 


0.16 


0.18 


0.22 


0.16 


0.13 


0.11 


1 










12:56:43.52 


+27:10:43.7 


5279 


52.4 


259.9 


2.7 


-2.22 


0.24 


0.118 


0.152 


1.00 


5.20 


-6.26 


-2.01 


4.84 


1.45 


3.73 


8.64 


1.48 


5.53 


0.155 


0.306 


4.96 


3.13 


2.81 


1.84 


0.85 


0.80 


0.19 


0.13 


0.005 


0.006 


0.10 


0.17 


0.20 


0.13 


0.26 


0.14 


0.21 


0.32 


0.14 


0.30 


0.003 


0.004 


0.14 


0.17 


0.20 


0.15 


0.11 


0.09 













12:55:41.30 


+27:15:02.7 


5886 


51.5 


258.3 


2.7 


-2.14 


-0.07 


0.110 


0.150 


1.02 


5.30 


-6.21 


-1.86 


5.00 


1.27 


3.46 


8.95 


1.49 


5.32 


0.160 


0.309 


4.91 


2.92 


2.78 


1.93 


0.80 


0.93 


0.20 


0.13 


0.005 


0.007 


0.10 


0.17 


0.20 


0.13 


0.26 


0.14 


0.22 


0.32 


0.14 


0.31 


0.003 


0.004 


0.14 


0.17 


0.21 


0.16 


0.12 


0.10 
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APPENDIX D: DATA TABLES 

The index and velocity dispersion measurements used in this work 
are tabulated in Table lDll 



